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Abstract We describe the current status of solar neutrino measurements and of the theory 
- both neutrino physics and solar astrophysics - employed in interpreting measurements. Im- 
portant recent developments include Super-Kamiokande's determination of the v — e elastic 
scattering rate for S B neutrinos to 3%; the latest SNO global analysis in which the inclusion of 



2 

low-energy data from SNO I and II significantly narrowed the range of allowed values for the 
neutrino mixing angle #12; Borexino results for both the 7 Be and pep neutrino fluxes, the first 
direct measurements constraining the rate of ppl and ppll burning in the Sun; global reanalyses 
of solar neutrino data that take into account new reactor results on #13; a new decadal evalua- 
tion of the nuclear physics of the pp chain and CNO cycle defining best values and uncertainties 
in the nuclear microphysics input to solar models; recognition of an emerging discrepancy be- 
tween two tests of solar metallicity, helioseismological mappings of the sound speed in the solar 
interior, and analyses of the metal photoabsorption lines based on our best current description 
of the Sun's photosphere; a new round of standard solar model calculations optimized to agree 
either with helioseismology or with the new photospheric analysis; and, motivated by the solar 
abundance problem, the development of nonstandard, accreting solar models, in order to inves- 
tigate possible consequences of the metal segregation that occurred in the proto-solar disk. We 
review this progress and describe how new experiments such as SNO+ could help us further 
exploit neutrinos as a unique probe of stellar interiors. 
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1 INTRODUCTION 



In 1958 Holmgren Sz Johnston (1958, 1959) found that the cross section for He + 



4 He 7 Be + 7 was about 1000 times larger than anticipated, so that in addition 
to the simplest 3 He + 3 He — > 4 He + 2p ppl termination of the pp chain (see Fig. 
[TJ, there might be significant branches to the ppll and ppIII cycles, and thus 
significant fluxes of 7 Be and 8 B solar neutrinos. Despite the uncertainties that 
existed in 1958 - the solar core temperature was poorly constrained by theory, 
and other nuclear physics important to the pp chain had not been resolved - 



both Cameron (1958) and Fowler (1958) pointed out that it might therefore be 



possible to detect solar neutrinos using a radiochemical method Ray Davis had 
developed at Brookhaven ( Davis Jr.|1955 ). While the endpoint of the main source 
of neutrinos from the ppl cycle, p+p— >d+e + +v e , is below the 811 keV threshold 



for ^ e + 37 Cl— > 37 Ar + e~, most 7 Be and 8 B neutrinos are sufficiently energetic to 
drive this reaction. In 1962 Fowler organized a team of young Caltech researchers 
- John Bahcall, Icko Iben, and Dick Sears - to begin the development of a solar 
model to more accurately predict the central temperature of the Sun and to 



estimate the rates of neutrino-producing reactions (Bahcall et al. 1963). The 



history of these early developments is summarized in several sources (Bahcall & 
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advances in the solar model and in the understanding of the nuclear physics of the 



pp chain and the 37 Cl(f e ,e~) Ar reaction, Davis Jr. (1964) and Bahcall (1964) 
concluded that a measurement of solar neutrinos would be possible, were Davis 
to mount a detector 100 times larger than that he built at Brookhaven, in a site 
sufficiently deep to reduce backgrounds from high-energy cosmic ray muons to an 



acceptable level. In April 1968 Davis Jr., Harmer k, Hoffman (1968) announced 
an upper bound on the solar neutrino capture rate for 37 C1 of 3 SNU (1 SNU = 
10~ 36 captures/ tar get atom/s), based on the initial running of a 100,000-gallon 
C2CI4 detector that the collaborators had constructed on the 4850- ft level of the 
Homestake gold mine, in Lead, South Dakota. 



This upper bound, nearly a factor of three below the rate predicted by the Bah- 



call, Bahcall & Shaviv (1968) standard solar model (SSM), began a controversy 



that took 30 years to resolve. Twenty of these years passed without independent 
confirmation of the Davis result: as the CI rate was a fraction of a count per 
day in 0.6 kilotons of organic liquid, other technologies with comparable sensi- 
tivity were not easily developed. Because the Davis experiment was sensitive 
to a flux of neutrinos that varies sharply with the solar core temperature (0( 8 B) 
~ Tq 2 ( Bahcall|l989 )), the result could be accommodated by a variety of possible 



changes in the SSM having the net effect of reducing Tq by ~ 5%. But as ad- 
ditional constraints on solar neutrino fluxes were established by the Kamiokande 
( Kamiokande Collaboration, K. S. Hirata et al.|1989| >, SAGE ( |SAGE Collaborfr 
tion, J. N. Abdurashitov et al.|1994| ), and GALLEX ( jGALLEX Collaboration, P 
Anselmann et aL||1994 ) collaborations, a more detailed pattern of fluxes emerged 
that was not easily reconciled with the SSM. In contrast, with the discovery of 
the MSW mechanism ( |Mikheyev fc Smirnov"||1985[ |1986| |Wolfenstein| 1 1 978a|b[ ) , 
it became apparent that neutrino oscillations augmented by matter effects could 
account for the observations, even for a small mixing angle. The conclusions of 



an Annual Reviews article from this period (Haxton 1995) captures the sense of 



excitement that with new experiments, a resolution of the solar neutrino problem 



was near. 



In 1998 —7- v T vacuum neutrino oscillations were discovered through the 
Super-Kamiokande collaboration's study of the zenith-angle dependence of at- 



mospheric neutrino fluxes (Super-Kamiokande Collaboration, Y. Fukuda et al. 



1998). While this result did not directly constrain the v e s produced by the Sun, 



the discovery was a game-changer, confirming a phenomenon originally suggested 



by Pontecorvo (1967) as a possible explanation for the solar neutrino problem. 



Finally, the Sudbury Neutrino Observatory (SNO) collaboration (SNO Collabo- 



ration, Q. R. Ahmad et al. 2001 2002) measured both the v e and heavy- flavor 



components of the solar neutrino flux arriving at Earth, utilizing three different 
detection channels with varying sensitivities to charge and neutral currents. The 
SNO collaboration measured the electron and heavy flavors components of the 8 B 
solar neutrino flux, found that the total flux of neutrinos (summed over flavors) 
is in good agreement with the SSM prediction, and determined flavor mixing pa- 



rameters that attributed the differential suppression of the pp, 7 Be, and 8 B fluxes 
deduced from previous experiments to the energy-dependent effects of matter on 
oscillations. 

This review summarizes the basic physics of solar neutrinos, the work that 
led to the discoveries noted above, and the impact of recent and ongoing solar 
neutrino experiments on astrophysics and weak interactions, including 



Completion of phase III of the Super-Kamiokande experiment (Super-Kamiokande 



Collaboration, K. Abe et al. 2011) and preliminary results from Super- 



Kamiokande IV's low-threshold running dSmy|2012| >; 



SNO's combined analysis of all three SNO phases (SNO Collaboration, B 



Aharmim et al. 2012) and the collaboration's low-energy analysis of the 



data from SNO I and II ( [SNO Collaboration, B. Aharmim et al.||2010| ); 



Borexino's achievement of a 5% measurement of the 7 Be flux, an initial re- 



sult for the pep flux, and a limit on the CN neutrino contribution (Borexino 



Collaboration, G. Bellini et al. 2011, 2012a); and 



Daya Bay, Reno, and Double Chooz measurements of #13, impacting global 



analyses of solar neutrino data (Daya Bay Collaboration, F. P. An et al. 



20T21 IDouble Chooz Collaboration, Y. Abe et al.||2012[ IRENO Collabora- 



tion, J. K. Ahn et al. 2012) 



In addition, a comprehensive survey of the nuclear physics of the pp chain and 
CNO cycle has been completed, yielding a new set of best values and uncertainties 
for the nuclear rates ( Adelberger et al.]|2011 ). The sound speed throughout most 
of the Sun's interior has been extracted from helioseismology to an accuracy 
~ 0.1%, providing a stringent check on SSM predictions. More sophisticated 
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3D models of the solar atmosphere have been developed, significantly improving 
the agreement between predicted and observed absorption line-shapes and the 
consistency of abundance determinations from different atomic and molecular 



lines (Asplund et al. 2009) - but also yielding a photospheric metal abundance 
~ 30% below ID values, leading to a conflict between SSMs employing the new 
abundances and solar parameters deduced from helioseismology. The SSM has 
been updated for the new nuclear reaction rates and alternative metallicities, 
and nonstandard models have been developed to explore accretion as a possible 
solution to the "solar abundance problem" ( Serenelli et al.|2009 Serenelli, Haxton 



& Peha-Garay 2011). 



For three decades solar neutrino physics was defined by an incomplete knowl- 
edge of the neutrino fluxes and shortcomings in our understanding of neutrino 
flavor physics. We are now starting a new period, where precise spectroscopy of 
the full spectrum of solar neutrinos is possible, and where a clearer understand- 
ing of weak interactions has been obtained from a combination of astrophysical, 
reactor, and accelerator experiments. On one hand, this returns us to the roots of 
solar neutrino physics: with weak interaction uncertainties removed, solar neutri- 
nos can be used to probe possible limitations in the SSM - such as uncertainties 
in the Sun's initial composition and the absence of ab initio treatments of mix- 
ing and other three-dimensional physics, including rotation and magnetic fields. 
On the other hand, the neutrinos coming from the Sun remain important to 
fundamental physics: the spectral shapes and fluxes of the various sources are 
known rather precisely, and low-energy neutrinos react with targets rather sim- 
ply, giving us confidence that we can interpret measurements. Thus this review 
also considers the continuing role solar neutrinos could play in further constrain- 
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ing the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) neutrino mass matrix and 
in probing matter effects and other environmental neutrino phenomena. 



2 THE SSM AND ITS NUCLEAR AND NEUTRINO PHYSICS 
2.1 The Standard Solar Model 

Solar models trace the evolution of the Sun from its beginning - when the collapse 
of the pre-solar gas cloud was halted by the turn-on of thermonuclear reactions 
- to today, thereby predicting contemporary solar properties such as the compo- 
sition, temperature, pressure, and sound-speed profiles and the neutrino fluxes. 
SSMs are based on four assumptions ( Bahcall||1989 ): 



The Sun burns in hydrostatic equilibrium, maintaining a local balance be- 
tween the gravitation force and pressure gradient. To implement this con- 
dition, an equation of state (EoS) is needed. As the hydrogen and helium in 
the Sun's core are nearly completely ionized, an ideal gas EoS with correc- 
tions for incomplete ionization of metals, radiation pressure, and screening 
is thought to provide a good approximation to the EoS ( Bahcall|1989 ). He- 



lioseismic inversions of solar p-mode frequencies have provided important 



tests of our understanding of the associated theory (Elliott & Kosovichev 



1998) 



The mechanisms for energy transport are radiation and convection. The in- 
ner portion of the Sun - 98% by mass or about 71% by radius - is radiative. 
In order to describe radiative transport, the opacity must be known as a 
function of temperature, density, and composition. In addition to elemen- 
tary processes such Thomson scattering off electrons and inverse bremm- 



strahlung off fully ionized hydrogen and helium, more complex processes 
such as bound-free scattering off metals are important contributors to the 
opacity in the Sun's central regions. Thus modern opacity tables like OPAL 



(Rogers & Nayfonov 2002) are based on detailed atomic input. Changes 
in opacity influence important helioseismic properties such as the sound 
speed and the location of the convective zone boundary. In the Sun's outer 
envelope, where the radiative gradient is larger, convection dominates the 
energy transport. SSM convection is modeled through mixing length the- 
ory, in which volume elements are transported radially over a characteristic 
distance determined empirically in the model, but typically on the order of 
the pressure scale height. 

The Sun produces its energy by fusing protons into 4 He, 

2e~ + 4p -> 4 He + 2v e + 26.73 MeV. (1) 

via the pp-chain (~ 99%) and CN cycle (~ 1%). The nuclear cross sec- 
tions are taken from experiment or, if that is impractical, from theory: the 
associated laboratory astrophysics is challenging because reaction rates are 
needed for energies well below the Coulomb barrier. Thus laboratory mea- 
surements are generally made at higher energies, with theory guiding the 
extrapolations to the solar Gamow window. 

Boundary conditions include the modern Sun's mass, age, radius R & , and 
luminosity L Q . The pre-solar composition is divided into hydrogen Xj n i, 
helium Yi n ;, and metals Zj n i, with Xi n i+Yi n i+Zi n i=l. Relative metal abun- 
dances can be determined from a combination of photospheric (determined 
from analyses of absorption lines) and meteoritic (for refractory elements) 
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abundances, and are generally assumed to have remained constant since 
the Sun formed. The photospheric abundances and the assumption of a 
homogeneous zero-age Sun then constrain the Sun's initial core composi- 
tion: one can equate the Sun's pre-solar core metallicity Zj n j to its present 
photospheric Zg, corrected for the effects of diffusion over the Sun's life- 
time. Finally Yi n j and the mixing length «mlt are determined interatively 
by demanding that L and Rq are produced after 4.6 Gy of burning. 

The resulting model is dynamic. The luminosity of the Sun increases by ~ 40% 
over the solar lifetime: Helium synthesis alters the mean molecular weight and 
opacity in the core, requiring a response in Tq- The ratio of ppl/ppll/pplll burn- 
ing changes, with the fraction of energy produced through the more temperature- 
sensitive ppll and ppIII branches increasing. The 8 B neutrino flux for the ppIII 
cycle has an exceedingly sharp growth ~ e t l T where r ~ 0.9 Gy. Composition 
gradients are created as the pp chain burns to equilibrium. An interesting exam- 
ple is the solar core 3 He abundance, which rises steeply with increasing radius, 
Xi oc T~ 6 , throughout the region where pp-chain equilibrium has been reached. 
The 3 He abundance gradient was proposed as a potential trigger for periodic 



mixing of the core in the "solar spoon" (Dilke & Gough 1972). Metals are rear- 



ranged: in the first 10 8 years of main-sequence burning most of the carbon in the 
Sun's central core is converted to nitrogen, building up the core abundance of 
14 N. Because they have a smaller charge-to-mass ratio than hydrogen, 4 He and 
metals slowly diffuse toward the core - another source of composition gradients 
that affect contemporary observables like helioseismology. 

Properties of two SSMs we will use in this review are listed in Table [TJ The 
models differ in the values assumed for the photospheric metallicity Zg, with the 
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GS98-SFII SSM being more metal rich than the AGSS09-SFII SSM. The table 
gives the model photospheric helium Y$ and metal Zg abundances, the radius of 
the convective zone Rcz, the mean deviation of the sound speed (5c/ c) from the 
helioseimic profile, the core helium and heavy element abundances Yc and Zq, 
and the Sun's pre-solar abundances Y; n i and Zj n i. 

2.2 The pp Chain and CN cycle 

Approximately 80% of the observed stars lie along a path in the Hertzsprung- 
Russell diagram characterized by energy generation through proton burning. The 
Sun provides a unique opportunity to test our understanding of such main- 
sequence stars, as we can compare model predictions to solar properties that 
are precisely known. This has inspired a great deal of laboratory work to reduce 
uncertainties in atomic opacities and nuclear cross sections - key SSM input pa- 
rameters - so that we can assess the reliability of the more fundamental solar 
physics and weak interactions aspects of the model. 

In lower-mass hydrogen-burning stars 4 He is synthesized primarily through the 
pp-chain nuclear reactions shown on the left side of Fig. [l| The rates of the ppl, 
ppll, and ppIII cycles comprising the pp chain can be determined from the fluxes 
of the pp/pep, Be, and 8 B neutrinos produced by those cycles. As the relative 
rates are very sensitive to Tc, the neutrino fluxes are a sensitive thermometer for 
the solar core, provided the associated nuclear physics is under control. 

Rates depend on the quantity (o"u)mb where v is the relative velocity between 
two colliding nuclei, a is the cross section, and ( )mb denotes an average over the 
Maxwell-Boltzmann relative velocity distribution in the solar plasma. The opti- 
mal energy for a solar reaction, called the Gamow peak, is determined from the 
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competition between a cross section that rises rapidly as the energy climbs the 
Coulomb barrier, and a relative-velocity distribution that declines rapidly on the 
Maxwell-Boltzmann tail. Two ppl-cycle reactions, d+p and 3 He+ 3 He, have been 
measured in their Gamow peaks. Data were obtained down to 2.5 and 16 keV, 
respectively, at LUNA, Gran Sasso's low-background facility for nuclear astro- 



physics (Bonetti et al. 1999). For other pp-chain reactions, direct measurements 
are not currently possible because of the severity of the Coulomb suppression. 
Instead, measurements must be made at higher energies, then extrapolated to 
solar energies using nuclear models to predict the cross section shape. 
Such extrapolations are usually performed by using the S-factor, 

a (E) = ^ exp [-2^(E)) , (2) 

which removes from the cross section the rapid energy dependence associated 
with the s-wave interaction of point nuclei. Here the Sommerfeld parameter 
r]{E) = Z\Zi a/v, with Z\ and Z2 the ion charges, the relative velocity v = 
y/2E/fi, [i is the reduced mass, and a is the fine structure constant (h = c = 
1). The remaining nuclear physics (including effects of finite nuclear size, higher 
partial waves, antisymmetrization, and any atomic screening effects not otherwise 
explicitly treated) is then isolated in S(E), the function used in extrapolations 
because of its gentler dependence on E. In solar applications S(E) is frequently 
approximated by its zero-energy value S(0) and corrections determined by its first 
and second derivatives, S'(0) and S"(0). 



The recent review by Adelberger et al. (2011 ) (Solar Fusion II) details the data 
and theory issues affecting our understanding of solar cross sections. Uncertain 
S-factors remain one of the key limitations in SSM neutrino flux predictions. 
Figure [2] gives the Solar Fusion II summaries for 3 He+ 3 He— )- 4 He+2p (left panel) 
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and 7 Be+p— )- 8 B+7 (right panel). While measurements for the first reaction cover 
the solar energies of interest, the screening environment of a terrestrial target 
(neutral atoms) differs from that at the center of the Sun (ionized 3 He). It is 
apparent from Fig. [2] that the theoretical screening correction is significant. 

The reaction 7 Be(p,7) 8 B (right panel of Fig. [2]) feeds the ppIII cycle that pro- 
duces the 8 B neutrinos measured by SNO and Super-Kamiokande. This reaction 
was considered the most uncertain in the pp chain when these detectors began 
operations, with only a single data set considered sufficiently reliable and well 



documented for inclusion in an S-factor determination (Adelberger et al. 1998). 
Four new, high-quality data sets were available for the Solar Fusion II evaluation, 
yielding Si7(0)=20.8 ± 0.7 (expt) ± 1.4(theor). The dominant error is now the 
theoretical extrapolation to the Gamow peak. 

The CN I cycle, illustrated on the right in Fig. [TJ is typically the dominant 
mode for hydrogen burning in massive main-sequence stars, where core temper- 
atures exceed those of the Sun. Unlike the pp chain, the CN cycle depends on 
pre-existing metals to catalyze a series of (p,7) and (p,a) reactions, leading in 
sum to Eq. 0. Thus the CN cycle is (in most settings) a secondary process, 
dependent on the star's metallicity. In the Sun the CN cycle converts C to N 
as it seeks equilibrium. Equilibrium has been reached only in the most central 
regions of the core, where T > 1.33 x 10 7 K. Outside this region, the bottleneck 
reaction 14 N(p,7) inhibits cycling. Thus, throughout most of the cooler regions of 
the core, very little CN-cycle burning takes place: pre-solar 12 C has already been 
converted to 14 N, but little 14 N is being consumed. Still further outward, where 
T < 10 7 K, the 12 C lifetime is comparable to the solar age. This is the region in 
the contemporary Sun where pre-solar 12 C is being burned. Thus CN-cycle neu- 
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trinos, produced in the /3 + decay of 15 O and 13 N, come from two distinct regions. 
Deep in the solar core equal numbers of 15 and 13 N neutrinos are produced in 
equilibrium burning, while in the distant cool outer regions of the core, only 13 N 
neutrinos are produced. 

2.3 Solar Neutrino Fluxes 

The main neutrino-producing reactions of the pp chain and CN cycle are summa- 
rized in Table [2] Four of the five (3 decay reactions produce neutrino spectra with 
allowed shapes and endpoints given by E™ ax . In the fifth, the decay of 8 B, the 
8 Be final state is a broad (~ 2 MeV) resonance. As the profile of this resonance 
is known, the associated small deviations from an allowed spectrum can be cal- 
culated. In addition to the main pp/pep, 7 Be, and 8 B neutrinos, a fourth source 
from a weak side-branch of the pp chain, the hep or 3 He+p neutrinos, is shown 
in Fig. [TJ These neutrinos are the most energetic produced by the Sun (E™ ax ~ 
18.77 MeV), and thus may be observable in SNO and Super-Kamiokande they 
populate energy bins above the 8 B neutrino endpoint. The two electron-capture 
(EC) reactions, p+e _ +p and Be+e~, produce line sources of neutrinos of energy 
E™ ax , broadened by thermal effects. There are two lines associated with Be EC, 
as ~ 10% of the captures populate the first excited state in 7 Li. 

The ppl, ppll, and ppIII contributions to solar energy generation can be de- 
termined from measurements of the pp/pep, 7 Be, and 8 B neutrino fluxes. As we 
will discuss in the next section, the Be and 8 B fluxes are now quite well known, 
while the first measurement of pep neutrinos was very recently announced. The 
"solar values" of Table [2] come from the Borexino Collaboration (private commu- 
nication) which updated an analysis by Bahcall Sz Pena-Garay| p003), combining 
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8 B, Be, and pep flux measurements (as available in March, 2011) with the solar 
luminosity constraint, to fix the principal pp-chain fluxes. That is, the sum of 
the rates for the ppl, ppll, and ppIII cycles, weighted by the energy each cycle 
deposits in the Sun, is fixed by the photon luminosity. Consequently the "solar 
values" are not strictly measured ones, but derived assuming a steady-state Sun. 

Table [2] also gives fluxes for two solar models, reflecting the metallicity un- 
certainties mentioned previously. The first model uses abundances for volatile 
elements that were obtained from an absorption line analysis in which the photo- 
sphere was treated in ID, yielding (Z/X) s = 0.0229 ( |Grevesse fc Sauval|19 98). As 
Solar Fusion II cross sections are used as well, this model is denoted GS98-SFII. 
The second, denoted AGSS09-SFII, takes abundances from a 3D photospheric 
model, yielding (Z/X) 5 =0.0178 ( |Asplund et aL1|2009D . The solar core is sensitive 
to metallicity, as free-bound/bound- free transitions in metals are an important 
contributor to the opacity. A low-metallicity Sun, as in model AGSS09-SFII, pro- 
duces a somewhat cooler core (by ~ 1%), and thus lower fluxes of temperature- 
sensitive neutrinos, such as those from 8 B (3 decay. 

The SSM fluxes for the CN I cycle /3 decays of 13 N and 15 O and the CNO II 
cycle decay of 17 F are also shown. Despite the minor role the CN cycle plays 
in solar energy generation, these fluxes are significant. The excess in the flux 
of 13 N neutrinos relative to 15 O neutrinos in Table 12] is a consequence of the 
out-of-equilibrium burning of the CN cycle discussed previously. 

The SSM uncertainties given in Table [T] are generated from the uncertainties 
assigned to approximately 20 model input parameters, denoted These in- 
clude the solar age, present-day luminosity, opacities, the diffusion constant, the 
S-factors for the pp chain and CN cycle, and the various metal abundances (key 
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elements such as C, N, O, Ne, Mg, Si, S, Ar, and Fe). The consequences of input 
SSM uncertainties on observables are typically parameterized through logarith- 
mic partial derivatives a(i,j), determined by calculating the SSM response to 
variations in the (3j. For example, in the case of the neutrino fluxes (pi, the 

ain[&M(o)] ^ 



can be found in the SSM updates of Peha-Garay & Serenelli (2008) and Serenelli 



( |2010p . Here ^(0) and /3j(0) denote the SSM best values. 

The partial derivatives define the power-law dependencies of neutrino fluxes 
with respect to the SSM best- value prediction ^(0), 

19 r a ia(i,j) 19 
3=1 lPj{ >j 3=1 

where 5(5 j = A(3j/(3j(0) is the fractional uncertainty of input parameter (5j with 
respect to its SSM best value. As this expression separates the impact of SSM 
parameter variations on <f>i into a solar piece - the infinitesimal SSM response 
described by a(i, j) - and a laboratory or theory piece - the estimated uncertainty 
5f3j of an input parameter, the effects of parameter variations can be explored 
without repeating SSM calculations. 

The solar abundance problem is characterized by large systematic differences 
in the SSM (3j for key abundances. Consequently the differences in the GS98- 
SFII and AGSS09-SFII SSM neutrino flux predictions of Table [2] exceed their 
respective internal statistical uncertainties, in several cases. The Table summa- 
rizes key helioseismic predictions of both models, with the low-Z AGSS09-SFII 
SSM predictions being in significantly poorer agreement with the data. 

The spectra of solar neutrinos emitted by the Sun is shown in Fig. [3j This 
familiar figure its somewhat idealized: it includes competing (3 decay and EC 
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branches for the p+p reaction, but omits the EC lines that accompany the other 
j3 decay reactions of Table [2] The EC branching ratio increases with increasing Z 
and decreasing Q- value. Thus, among the omitted cases, EC is most significant for 



the CNO cycle reactions, where rate estimates were made by Stonehill, Formaggio 



& Robertson (2004). The EC lines shown in the figure are in fact thermally 



broadened, as they occur in a plasma. Finally, at energies < 10 keV below the 
scale of Fig. |3j there is a contribution from neutrino pairs of all flavors produced 
thermally in the solar core ( Haxton Sz Lin|2000 ) : while the flux of these neutrinos 
is modest, the peak flux density of ~ 10 9 /cm 2 /s/MeV exceeds that of all solar (3 
decay sources, except for the pp neutrinos. 

In the following sections, measurements of the various sources can provide 
unique information on the structure and composition of the Sun, and on the 
properties of neutrinos, including how those properties depend on the surrounding 
environment. 



3 EXPERIMENTS: NEUTRINOS AND HELIOSEISMOLOGY 

The SSM is a model of the Sun's interior structure, rather than the more compli- 
cated behavior of the convective envelope. The two main tools by which we can 
probe the solar interior and thus test the validity of the SSM are neutrino spec- 
troscopy and helioseismology. Neutrino fluxes are sensitive to core temperature 
- a property that responds to changes in the opacity and composition - and to 
nuclear cross sections. Helioseismic maps of the sound speed throughout much of 
the solar interior have achieved accuracy of a few parts in 1000, constraining the 
solar density and pressure profiles, and determining rather precisely the boundary 
between the Sun's radiative core and convective envelope. 
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3.1 Neutrino Spectroscopy: Early Experiments 

3.1.1 The Chlorine Experiment Radiochemical detection of neutrinos 
by 37 Cl(V e ,e _ ) 37 Ar was suggested by Pontecorvo (1946) and explored in more 



detail by |Alvarez ( 1949 ), who was interested in the possibility of a reactor neutrino 



experiment to detect a Majorana neutrino. Ray Davis's efforts to develop a 
practical detector began with his Brookhaven experiment ( Davis Jr.|1955 ), which 



used a 1000-gallon tank of perchloroethylene (C2CI4) placed 19 feet underground. 
This yielded an upper bound on the rate for solar neutrino reactions of ~ 40,000 



SNU. Subsequent developments are described by Lande (2010). Construction 
began on the Homestake detector in 1965, with first results announced in 1968, 
and with measurements continuing until 2002, when the Homestake Mine closed. 
The final result 

{&(()) = 2.56 db 0.16 ± 0.16 SNU (5) 

is about a factor of three below historical SSM best values. (The SSM GS98-SFII 
rate is 8.00 ± 0.97 SNU.) 

The experiment exploited fortuitous properties of 37 Ar in achieving nearly 
single-atom counting. The average solar neutrino reaction rate in the tank 
was 0.48 counts/day, above a background dominated by cosmogenics of 0.09 
counts/day. As a noble gas that does not interact chemically, argon can be 
extracted with high efficiency (> 95%) from large volumes of organic liquid. The 
37 Ar half life of ~ 35 days allowed tank concentrations to build up over a satu- 
ration time of ~ two months, yet also made 37 Ar counting via electron capture 
feasible. As the decay populates the ground state of 37 C1, the signal is the 2.82 
keV Auger electron produced as the electrons in 37 CI adjust to fill the K-shell va- 
cancy. Davis developed miniaturized gas proportional counters for this counting. 
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Taking into account detector efficiencies and losses due to 37 Ar decaying in the 
tank, the number of Ar atoms counted per year was ~ 25. 

The chlorine experiment was primarily sensitive to the temperature-dependent 
neutrino fluxes produced in the ppIII and ppll cycles ( 8 B ~ 75%, 7 Be ~ 16%). 
For this reason the source of the "solar neutrino problem" was not immediately 
apparent. Many candidate explanations were offered over the next 30 years, with 
many of these proposing changes in the SSM to produce a somewhat cooler core. 
3.1.2 Kamiokande II/III Confirmation of the 37 C1 results came 21 years 
later, from a detector originally designed for proton-decay studies. The Kamiokande 
I proton decay detector was upgraded in the early 1980s to Kamiokande II/III, 
a 3.0-kiloton imaging water Cherenkov detector capable of detecting solar and 
other low-energy neutrinos. The neutrino signal is the Cherenkov light emitted 
by recoiling electrons after elastic scattering (ES), v x + — >■ v' x + e~, which is 
sensitive to both electron and heavy-flavor neutrinos, though with the differential 
sensitivity cr(i/ e )/cr(i/ /1 ) ~ 6. For incident neutrino energies 3> m e c 2 , the electron 
is scattered in the forward direction. Thus, by correlating event directions with 
the position of the Sun, one can cut away a large background uncorrelated with 
solar position, to reveal solar neutrino events in a forward cone. 

The inner 2.14 kilotons of the detector was viewed by 948 20-inch Hamamatsu 
photomultiplier (PMT) detectors, providing ~ 20% coverage. The innermost 0.68 
kilotons of the detector served as the fiducial volume for solar neutrino detection. 
Kamiokande II operated with a ~ 9 MeV threshold, which was later reduced to 
7.5 and 7.0 MeV in Kamiokande III. 

The improvements made in Kamiokande II to enable solar neutrino detection 
included water purification to reduce low-energy backgrounds associated with 
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radon and uranium as well as electronics upgrades to improve timing, vital for the 
reconstruction of the interaction vertices and directions of low-energy electrons, 
and thus in more cleanly defining a fiducial volume for solar neutrino events. 
After water-sealing the cavity holding the main detector, the outer portion of 
the detector was instrumented with 123 PMTs to serve as a muon anti-counter, 
and additional water was added to shield against 7s from the surrounding rock. 
Kamiokande III included improvements in the electronics, water purification sys- 
tem, event reconstruction and selection tools, and Monte Carlo detector simula- 
tions software. 

The first production run of Kamiokande II began in January 1987. The de- 
tection of 8 B solar neutrinos based on 450 days of data extending through May 



1988 was announced by Kamiokande Collaboration, K. S. Hirata et al. (1989). 
The measured flux of 8 B neutrinos with energies above 9.3 MeV was found to 
be 0.46 ± 0.13(stat) ± 0.08(syst) of the SSM value, confirming the deficit seen 
by Davis and collaborators. Kamioka II/III ran until February 1995, collecting 



2079 days of data. The combined analysis of all data sets yielded (Kamiokande 



Collaboration, Y. Fukuda et al. 1996) 



</>( 8 B) = (2.80 ± 0.19(stat) ± 0.33(sys)) x 10 b /Ws, (6) 

or 50% (61%) of the GS98-SFII (AGSS09-SFII) SSM result. 

The Kamiokande II/III detector was the first to record solar neutrinos event 
by event, establish their solar origin through correlation with the direction to 
the Sun, and provide direct information on the 8 B spectrum through the recoil 
electron spectrum from ES. 

3.1.3 The Gallium Experiments Two radiochemical gallium experiments 
exploiting the reaction 71 Ga(i/ e ,e~) 71 Ge, SAGE and GALLEX/GNO, began so- 



21 

lar neutrino measurements in January 1990 and May 1991, respectively. SAGE, 
which continues to operate, uses a target of 50 tons of Ga metal, heated to 30° C 
so that the metal remains molten, and has reported results for 168 extractions 
through December 2007. The experiment is located in the Baksan Neutrino Ob- 
servatory, under Mt. Andyrchi in the Caucasus. GALLEX, which used 30 tons 
of Ga in the form of a GaCi3 solution, operated through 1997, and its successor 
GNO continued through 2003. GALLEX and GNO were mounted in the Gran 
Sasso National Laboratory, near L'Aquila, Italy. 



Gallium, first proposed as a solar neutrino detector by Kuzmin (1966), has 
a low threshold for solar neutrino absorption (233 keV) and a strong Gamow- 
Teller transition to the ground state of 71 Ge. This leads to a large cross section 
for absorbing the low-energy pp neutrinos. As 71 Ge has a half life of 11.43 days, 
a radiochemical experiment analogous to that done for chlorine was proposed, 
though the chemistry of Ge recovery is considerably more complicated than that 
for Ar. Because of its pp neutrino sensitivity, the Ga experiment has a minimum 
astronomical counting rate of 79 SNU, assuming only a steady-state Sun and 
standard-model (SM) weak interaction physics. That is, any combination of ppl, 
ppll, and ppIII burning consistent with the Sun's observed luminosity will lead 



to a solar neutrino capture rate at or above this value (Bahcall 1989). Thus the 
gallium experiment had the potential to yield a result that would require a "new 
physics" solution to the solar neutrino problem. 

In 1974 Davis and collaborators began work on the chemistry of Ge recovery for 
both GaCi3 solution and Ga metal, conducting a 1.3-ton pilot experiment using 
GaCbj in 1980-82 to demonstrate the procedures later used by GALLEX. The 
method recovers Ge as GeCL: by bubbling nitrogen through the solution, then 
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scrubbing the gas. The Ge can be further concentrated and purified, converted 
into GeH4, then counted in miniaturized gas proportional counters similar to 
those used in the chlorine experiment. 

In the SAGE experiment with room-temperature liquid metal, the 71 Ge is 
separated by mixing into the gallium a solution of hydrogen peroxide and dilute 
hydrochloric acid, which produces an emulsion, with the germanium migrating 
to the surface of the emulsion as droplets, where it can be oxidized and dissolved 
by hydrochloric acid. The Ge is extracted as GeCU, purified and concentrated, 
synthesized into GeH4, then counted as in the GALLEX experiment. In both 
GALLEX and SAGE, the overall efficiency of the chemical procedures can be 
determined by introducing Ge carrier. 

A unique aspect of the gallium experiments was the neutrino source exper- 
iments done to check overall experimental procedures - chemical extraction, 
counting, and analysis techniques. The calibrations also checked the capture 
cross section, as two excited-state transitions not constrained by the 71 Ge life- 
time contribute to 7 Be neutrino capture. Two GALLEX calibrations and the 
first SAGE calibration were done with 51 Cr sources, while the second SAGE cal- 
ibration used an 37 Ar source. Source intensities were ~ 0.5 MCi. The weighted 



average (SAGE Collaboration, J. N. Abdurashitov et al. 2009) for the four cali- 
brations, expressed as the ratio R of measured 71 Ge to that expected due to the 
source strength, is R = 0.87 ± 0.05 (la). The discrepancy, which exceeds two 
standard deviations, has attracted some attention due to other short-baseline 



neutrino anomalies (Abazajian et al. 2012, Gavrin et al. 2011). 

SAGE began taking data in December 1989. The capture rate limit obtained 
from five extractions with 30 tons of gallium, < 79 SNU (90% c.l.) ( |SAGE Col- 
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laboration, A. I. Abazov et al. 1991), coincided with the minimum astronomical 



value. The most recent SAGE combined analysis for all 168 extractions yielded 



(SAGE Collaboration, J. N. Abdurashitov et al. 2009) 



(visage = 65.4tl;i(stat)^;°(syst) SNU 



+2.6/ 



(7) 



or approximately half the un-oscillated SSM best value. GALLEX began taking 
data in May 1991 and announced first results a year later, a counting rate based 



on 14 extractions of 83 ± 19 (stat) ± 8 (syst) SNU (GALLEX Collaboration, 



P. Anselmann et al. 1992). GALLEX completed four campaigns, I through IV, 



running until 1997. From 65 extractions and 1594 days of data, and including 



updates due to new pulse-shape analysis methods, GALLEX found (GALLEX 
Collaboration, W. Hampel et al.||"l999[ |Kaether et"al]|2010 ) 



(ct<P)gallex i-iv = 73.li^(stat)^;{(syst) SNU. 



n+3.7/ 



(8) 



A number of improvements in Ge extraction procedures, electronics, counter effi- 
ciency calibrations, and radon event characterization were incorporated into the 
follow-up experiment GNO. The experiment accumulated 1687 days of running 
between May 1998 and April 2003. The counting rate from the 58 extractions is 



(^)gno = 62.9i^(stat)+^(syst) SNU. 



The weighted average of SAGE, GALLEX, and GNO results is 



(9) 



{&<i>) SAGE+GALLEX+GNO = 66.1 ± 3.1 SNU, 



(10) 



with all uncertainties combined in quadrature (SAGE Collaboration, J. N. Ab- 



durashitov et al. 2009). The SSM GS98-SFII rate is 126.6 ± 4.2 SNU. 



3.1.4 Hints OF New Physics In Fig. [4] results of the early experiments are 
compared to the predictions of the contemporary GS98-SFII SSM. Not only are 
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the results in disagreement with the SSM, but the pattern of discrepancies is not 
easily reproduced even if one entertains the possibility of substantial variations in 
that model. By the early 1990s several analyses ( Hata &: Langacker||1994 Parke 



1995; White, Krauss & Gates 1993) had pointed to apparent contradictions in 
the pattern of fluxes with respect to SSM predictions, 

0(pp) ~ O.90 SSM (pp) </>( 7 Be)~0 ct>{ s B) ~ 0.4</> SSM ( 8 B). (11) 

Now variations in the SSM affect the neutrino fluxes principally through their 
impact on the core temperature Tq- As 

<K 8 B) 



T c \ (12) 



<Kpp) 

the observation that </>( 8 B)/</>(pp) ~ 0.4(/> SSM ( 8 B)/(/> SSM (pp) would seem to re- 
quire a cooler solar core, Tc ~ 0.95 T C SM . On the other hand, as 



0( 7 Be) 12 
0(«B) ~ c ' 



Tc 12 , (13) 



the observation that ^( 7 Be)/c/>( 8 B) << <^ SSM ( 7 Be)/</> SSM ( 8 B) would seem to re- 
quire a hotter core, Tc > T,S SM , a contradiction. Model- independent analyses 
assuming undistorted neutrino spectra and a steady-state Sun (so that neutrino 



fluxes are constrained by the Sun's luminosity) were done by|Hata, Bludman fc 



Langacker (1994) and Heeger Sz Robertson (1996). Their calculations showed 



that the probability of solutions without new particle physics were in a range 



of ~ 2-3%. Heeger & Robertson (1996) further argued that if the luminosity 



constraint were relaxed, this probability would still be limited to ~ 4%. The 
likelihood of a new-physics solution to the solar neutrino problem was high. 

The conclusion - that the solar neutrino problem might have its origin outside 
of astrophysics - was additionally supported by a growing body of evidence from 
helioseismology that validated SSM descriptions of the Sun's interior structure. 
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3.2 Helioseismology 

Measurements and analysis of Doppler shifts of photospheric absorption lines 
show that the Sun's surface oscillates with amplitudes ~ 30m and velocities ~ 
0.1 m/s, reflecting a variety of interior modes ( Gizon &: Birch|2005 ) . Turbulence 



within the Sun's convective zone acts as a random driver of sound waves prop- 
agating through the gas. Specific frequencies are enhanced as standing waves, 
normal eigenmodes that reflect details of solar structure. Here we summarize 



the basics of solar oscillations, referring readers to Chaplin & Miglio (2013) for a 
more detailed discussion. 

The SSM is characterized by quasi-static pressure p(r), density p(r), tem- 
perature T(r), entropy s(r), gravitational potential 4>{r), and nuclear energy 
generation e(r) profiles that are functions of the radial coordinate r. One can 
perturb the SSM by introducing small displacements 5f and associated velocities 



v(r) = d5r/dt, then seek small-amplitude normal-mode solutions ( Christensen- 
Dalsgaard|[2002l > 



p(r, t) = p (r) + p'(r, t) p'(f, t) ~ p'(r) Y lm (d, 0)e^ (14) 

that might account for observed solar surface oscillations. 

Solar oscillations can be treated in the adiabatic approximation because the 
timescale for heat exchange is much longer than the oscillation periods of interest. 
From the adiabatic index T\ describing the power-law dependence of the pressure 
on the density and the associated sound speed c(r), 

d log p{r) \ 1 2 



r ' S laraJ, P( r ) = jTP(r)c(r), (15) 
one can define an auxiliary field ^(r) = c 2 y/ p(r) In the Cowling (neglecting 



perturbations to the gravitational field) and adiabatic approximations (Duebner 
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fc Gough|[l984l ) 

d 2 ^i{r) 1 



dr 2 c 2 



(16) 

where propagating (evanescent) solutions exist for ui 2 s > (<0). This eigenvalue 
problem is governed by the buoyancy, or Brunt- Vaisala, frequency N(r), 

2 _ Gm(r) ( 1 dlog p(r) (flog p(r) \ 
[r) r [V! dr dr J' 1 ' 

which turns negative in the convective zone but is positive and roughly constant 
in the radiative interior; the Lamb frequency, 

which diverges for r — > if I > 0; and the acoustic cutoff frequency, which 
depends on the density scale height H(r) and sound speed, 



^w-^-'f whcrc ms -{m* ] ' (19) 



and determines the outer turning point where U) ~ uj co . Eigensolutions of Eq. ( 16 ) 
can be found for discrete frequencies {uj n i}, where n is the radial order: there is 
no dependence on m because all azimuthal modes for fixed n, I are degenerate by 



spherical symmetry. The assumptions leading to Eq. ( 16) can be justified except 



when n and I are small, or when I <C n ( Duebner Sz Gough"]|1984 ). 



As w > lj co everywhere except near the surface, the solar regions supporting 
propagating solutions are determined by 



4~«"-^(l-£')>0. (20) 



Two different families of solutions exist. The g-mode family is determined by the 
conditions uj 2 -C N 2 and to 2 <C Sf. Consequently g-mode propagation is confined 
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to the solar radiative interior. The second family, the acoustic oscillations or 
p-modes, are the modes that have been observed in the Sun. If uj 2 3> N 2 , then 



Eq. (20) and the requirement uj 2 s > define the inner turning-point radius 



c(r) 



'"turning 



UJ 



VW+lj. (21) 



Qualitatively it is clear that the dependence of the eigenfrequencies on I can 
provide localized sensitivity to c(r), with modes of low I penetrating more deeply 
into the solar interior. Because the eigenfrequencies depend on c(r), the p-mode 
observations constraint the solar pressure and density profiles. 

Similar radial sensitivity is found for the g-modes. The condition uj <C N(r) 



allows propagation in the deep interior, as Eq. (20) guarantees that uj 2 s > 
for r sufficiently small. While in principle this suggests sensitivity to c(r) in the 
solar core, g-modes are damped in the convective envelope, making observation 



difficult. No undisputed detection exists to date (Appourchaux et al. 2010). 

The significant effort invested in helioseismological measurements and analysis 
has yielded a rather precise map of c(r) over the outer 90% of the Sun by radius. 
The solar profile used in Fig. [5] was obtained by Basu et al. ( 2009 ) from an 



analysis that included 4752 days of BiSON data (http://bison.ph.bham.ac.uk). 
The comparison SSMs are AGSS09-SFII ((Z/X)<? = 0.0178) and GS98-SFII 
((Z/X) s = 0.0229). 

GS98-SFII is representative of models that were in use in the 1990s: the gen- 
erally good agreement with the solar c(r) (~ 0.2% apart from a narrow region 
just below the convective boundary) was taken as strong support for the SSM, 
helping to reinforce the conclusion that the solar neutrino problem might have a 
particle physics origin. Helioseismic data forced improvements in the SSM, such 



as the inclusion of helium and heavy-element diffusion (Bahcall & Pinnsoneault 
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1995, Bahcall & Pinsonneault 1992). The suggestion from early solar neutrino 
experiments and helioseismology that new particle physics could be the source 
of the solar neutrino problem provided additional motivation for a new genera- 
tion of sophisticated experiments with high statistics and sensitivity to neutrino 
flavor, spectral distortions, and day-night differences, described below. 

3.3 Super-Kamiokande 

Super-Kamiokande, the successor to the Kamiokande detector, is a 50-kton cylin- 
drical water Cherenkov detector located in the Kamioka Mine at a depth of ~ 
2.03 km of water (flat site equivalent). The inner 32-ktons of water is viewed 
by ~ 11,100 20" photomultipliers (40% coverage), with 22.5-ktons serving as 
the fiducial volume for detecting solar neutrinos. The detector has operated at 
(total) electron energy thresholds for solar neutrinos ranging from 7.0 to the 
present 4.0 MeV, so that detection is limited to 8 B and hep neutrinos. While 
in ES the energy of the incident neutrino is shared between the scattered elec- 
tron and outgoing neutrino, electron detection provides some sensitivity to the 
initial neutrino spectrum and thus to distortions associated with neutrino oscilla- 
tions ( Super-Kamiokande Collaboration||2012 Super-Kamiokande Collaboration 



S. Fukuda et al. 2003). The electron energy resolution at 10 MeV is ~ 16%. 

The detector began operations in 1996, progressing from phase I to the cur- 
rent phase IV. Super-Kamiokande I recorded neutrino events for approximately 
five years, determining an 8 B neutrino flux of </>( 8 B) = (2.35 ± 0.02 (stat) ± 
0.08 (syst)) x 10 6 /cm 2 /s from events recorded above 5 MeV, assuming an undis- 



torted spectrum (Super-Kamiokande Collaboration, J. Hosaka et al. 2006). The 



measured rate variation of ~ 7% over the year is consistent in magnitude and 
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phase with the effects of the Earth's orbital eccentricity, 1.7%. No evidence was 
found for spectral distortions or day- night differences, two signatures of neutrino 
oscillations in matter. 

The detector was drained following phase I for repairs and maintenance. Dur- 
ing refilling the implosion of a phototube led to a blast wave that destroyed 
most of the lower part of the detector. Super-Kamiokande II was subsequently 
rebuilt with the remaining phototubes enclosed in blast shields. Despite the 
reduced phototube coverage of 19% and resulting higher threshold of 7 MeV, 
Super-Kamiokande II ran successfully as a solar neutrino detector for three years, 
beginning in late 2002. The deduced rate, </>( 8 B) = (2.38 ± 0.05 (stat) j^Jjj (sys)) 
x 10 6 /cm 2 /s, is consistent with Super-Kamiokande I results. No spectral distor- 
tion was detected, and the day-night difference was again consistent with zero at 
la ( Super-Kamiokande Collaboration, J. P. Cravens et al.||2008 ). 

Super-Kamiokande III ( Super-Kamiokande Collaboration, K. Abe et al.||2011 ) 
collected nearly two years of data between October 2006 and August 2008, operat- 
ing with a fully restored set of 11129 PMTs equipped with blast shields, providing 
40% phototube coverage. Improvements made to the water purification system, 
event reconstruction and selection tools, and the Monte Carlo detector simulation 
software resulted in a reduced systematic uncertainty of ±2.1%. The observed 
event rate for electrons between 5.0 and 20 MeV is equivalent to an unoscillated 
8 B neutrino flux of (2.39 ± 0.04 (stat) ±0.05 (sys)) x 10 6 /cm 2 /s flSmy|2012| >. No 
significant spectral distortion was observed. 

Preliminary results from 1069 days of running for Super-Kamiokande IV were 



reported at Neutrino 2012 (Smy 2012). This latest phase of Super-Kamiokande 
includes new electronics, an improved Monte Carlo model of the trigger efficiency, 
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higher efficiency due to relaxed cuts against radioactivity backgrounds, and a 
lower threshold of 4 MeV (total energy). 



3.4 Sudbury Neutrino Observatory 

A second remarkable detector, the Sudbury Neutrino Observatory (SNO), was 
constructed two kilometers below ground, within the INCO Creighton nickel 
mine, Ontario, Canada ( Jelley, McDonald &; Robertson|2009 ) . A kiloton of heavy 



water was contained in a 12m-diameter spherical acrylic vessel. A surrounding 
array of 9500 20-cm PMTs viewed the inner volume, providing ~ 56% coverage. 
Seven kilotons of light water provided a 5m buffer between the central detector 
and the surrounding rock walls. 



SNO was proposed by Chen (1985), who recognized the advantages of the 
multiple detection channels that could be introduced by replacing the hydrogen 
in an ordinary water Cherenkov detector with deuterium. The flux of higher 
energy solar electron neutrinos can be probed with the charged-current (CC) 
reaction 

v e + d^p + p + e~, (22) 

with detection of the produced electron. As the Gamow- Teller strength is con- 
centrated near the 1.44 MeV breakup threshold for deuterium, the electrons carry 
off most of the energy, and thus provide significant information on the incident 
neutrino spectrum. A second channel, the neutral-current (NC) reaction 

v x + d — > v' x + n + p, (23) 

is independent of the neutrino flavor x, counting all neutrinos above the 2.22 
MeV breakup threshold. As the only detectable signal of the reaction is the pro- 
duced neutron, this channel placed very stringent constraints on the radioactive 
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cleanliness of the detector. The third channel is the ES reaction of conventional 
water detectors, 

v x + e~ -»• v' x + e~'. (24) 

Operations were carried out in three phases. SNO I operated with pure heavy 
water. The NC-channel neutrons can capture on deuterium, producing 6.25- 
MeV gamma rays that Compton scatter off electrons, producing light for recoils 
above the Cherenkov threshold. SNO I operations covered 306.4 live days from 
November 1999 through May 2001. Two analyses were performed based on the 
assumption of an undistorted 8 B spectrum, using electron kinetic energy thresh- 
olds of 6.75 and 5 MeV, respectively. The second analysis thus included NC 
events. 

In SNO II two tons of purified NaCl were dissolved in the water, so that 
35 Cl(n,7) would become the dominant neutron sink. This reaction increases the 
capture rate and energy release. Data were accumulated for 391.4 live days from 
July 2001 through August 2003. Detector calibrations completed in SNO I were 
repeated and extended in SNO II, including new checks involving the introduction 
of beta-gamma sources that could lead to photo-disintegration of deuterium and 
the use of a 252 Cf neutron source to determine the neutron detection efficiency. 
The analysis was performed for a kinetic energy threshold of 5.5 MeV and treated 
the first 254.2 live days of data as blind. In addition, the 8 B spectrum shape was 
not assumed, but rather extracted from the analysis, using 0.5 MeV bins from 
5.5 to 13.5 MeV, plus an additional bin for events between 13.5 and 20 MeV. 

In the first two phases of SNO the CC, ES, and NC rates were determined by a 
statistical analysis that decomposed the common signal, the Cherenkov light, into 
the three contributing components. The analysis exploited distinguishing angular 
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correlations with respect to the Sun and energy differences in the CC-, ES-, and 
NC-associated light. In SNO III the separation of the NC and CC/ES signals 
was done by direct counting of NC neutrons. The salt introduced in SNO II was 
removed by reverse osmosis and ion exchange, and a month of data was taken to 
confirm that the detector had been restored to the operating conditions of SNO 
I. Then an array of the specially designed 3 He- and CF4-filled gas proportional 
counters was installed for neutron detection by 3 He(n,p) 3 H. This neutral-current 
detection (NCD) array consisted of 40 strings of proportional counters, ranging 
in length from 9 to 11 meters, that were anchored to the inner surface of the 
acrylic vessel, forming a lattice on a one-meter grid. 

Between November 2004 and November 2006 385.17 live days of SNO III data 
were taken. Extensive calibrations of both the NCD and PMT arrays were made, 
utilizing various neutron and gamma-ray sources, in order to calibrate the effec- 
tiveness of the neutron detection and the impact of array installation on detector 
behavior. The array was exploited to characterize neutron backgrounds within 
the detector, including the distribution and isotopic composition of background 
sources. During solar neutrino running, data were culled to eliminate strings 
that exhibited mechanical or electrical faults, or runs (operational periods of at 
least 30 minutes) when any array abnormalities were observed. A blind analysis 
of the remaining data was then performed. The neutrino spectrum was again 
determined from the CC and ES data, not assumed. 

The SNO I/II and SNO III results are in generally good agreement, and both 
separately and in combination established 

1. A total flux of active neutrinos from 8 B decay of </>nc (^active) = (5.25 ± 
0.16(stat)l '^3(syst)) x 10 6 /cm 2 /s, in good agreement with SSM predictions, 
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and 0cc(^e) ~ 0.34 </> NC (^active); and 

2. The absence of statistically significant day-night effects or spectral distor- 
tions in the portion of the 8 B neutrino spectrum above ~ 5 MeV; 

3.5 Borexino 



The Borexino experiment (Borexino Collaboration, G. Alimonti et al. 2009), lo- 



cated in the Gran Sasso Laboratory at an effective depth of about 3.0 km.w.e., 
is the first to measure low-energy (< 1 MeV) solar neutrino events in real time. 
The detector is housed within a 16.9m domed tank containing an outer layer 
of ultrapure water that provides shielding against external neutrons and gamma 
rays. At the inner edge of the water a stainless steel sphere serves as a support 
structure for an array of photomultiplier tubes that view both the inner detector 
and the outer water shield, so that the Cherenkov light emitted by muons pass- 
ing through the water can be used to veto those events. Within the steel sphere 
there are two regions, separated by thin nylon vessels, containing high-purity 
buffer liquid, within which is sequestered a central volume of 278 tons of organic 
scintillator. The fiducial volume consists of ~ 100 tons of the liquid scintillator 
at the very center of the detector. Scintillation light produced by recoil electrons 
after ES events is the solar neutrino signal. The 862 keV 7 Be neutrinos produce 
a recoil electron spectrum with a distinctive cut-off edge at 665 keV. 

The Borexino Collaboration reported results in 2008 and 2011 constraining the 



fluxes of three low-energy solar neutrino branches (Borexino Collaboration, G. 



Bellini et al. 2011 2012a): 



1. A 7 Be solar rate equivalent to an unoscillated flux of (3.10±0.15) x 10 9 /cm 2 s, 
or about 62% of the GS98-SFII SSM central value: 
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2. A ES rate for 8 B neutrinos, based on a integration above 3 MeV, corre- 
sponding to an equivalent flux of ES ( 8 B) = (2.4 ± 0.4 ± 0.1) x 10 6 /cm 2 s 



(Borexino Collaboration, G. Bellini et al. 2010), less precise than but in 



good agreement with SNO and Super-Kamiokande results. [A similar re- 
sult has been obtained by the KamLAND collaboration, <^> ES ( 8 B) = (2.77 ± 
0.26 ±0.32) x 10 6 /cm 2 s from events above their 5.5 MeV analysis threshold 
QKamLAND Collaboration, S. Abe etaIp0lTl >.] 



3. The first direct, exclusive determination of the pep flux, (1.6 ± 0.3) x 
10 8 /cm 2 s (95% c.l.); and 

4. A limit on the CNO neutrino flux, </>cno < 7.7 x 10 8 /cm 2 s at 95% c.l. 

The Borexino 7 Be measurement places an important constraint on matter effects 
in neutrino oscillations, as this line lies in a region dominated by vacuum os- 
cillations, while the Super-Kamiokande and SNO measurements are done in the 
matter-dominated region. 

4 NEW NEUTRINO PROPERTIES 

The results just described have been addressed in global analyses that extract 
from the experiments constraints on neutrino and solar properties. Before de- 
scribing such analyses, we discuss some of the associated weak interactions is- 
sues. As previously noted, by 1994 Kamiokande II /III had finished operations, 
confirming the neutrino deficit that Davis, Harmer, and Hoffman had first discov- 
ered 26 years earlier, and the SAGE and GALLEX experiments had converged 
on a counting rate very close to the minimum astronomical value of 78 SNU. The 
pattern of pp, 7 Be, and 8 B fluxes that emerged from analyses of the three early 



experiments (see Eq. ( 11 )) was inconsistent with possible SSM variations altering 
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Tc (see figures by Castellani et al. (1994) and Hata, Bludman Langacker (1994) 



included in Haxton (1995)) and improbable in model- independent analyses that 



assumed only undistorted neutrino spectra (Hata, Bludman & Langacker 1994 



Heeger &; Robertson|[X9 96 ) , The agreement between the SSM sound speed profile 



and that deduced from helioseismology also made it more difficult to motivate 
SSM changes. 

A variety of new particle physics solutions to the solar neutrino problem had 
been suggested over the years, including vacuum and matter-enhanced neutrino 



oscillations, neutrino decay (Bahcall, Cabibbo Ik Yahil 1972), and weakly in- 
teracting massive particles (WIMPs) that might be bound in the Sun and con- 



sequently contribute to energy transport (Faulkner & Gilliland 1985, Spergel 



& Press 1985). In addition to the standard MSW scenario, other oscillation ef- 



fects in matter were explored, including spin-flavor resonances driven by neutrino 
magnetic moments ( Akhmedov|| 1988" Lim &: Marciano]|1988 ), parametric density 
fluctuations ( Krastev &: Smirnov| [l989, Sc hafer Koonin"||1987 ), contributions to 
the MSW potential from currents (Haxton & Zhang 1991), and depolarization in 
the stochastic fluctuation limit ( Loreti Balantekin||1994 ). 

Of these and other possibilities, the MSW mechanism drew the most attention 
because of its minimal requirements, neutrino masses and a vacuum mixing angle 
9 V > 10 -4 . While neutrinos are massless in the SM, and consequently cannot 
oscillate, nonzero masses arise in most extensions of the model. Small weak 
interaction mixing angles were already familiar from the analogous quark mixing 
matrix. 
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4.1 Oscillation Basics: The Vacuum Case 



The current laboratory (tritium (3 decay) limit on the v e mass is 2.3 eV, though an 



effort is underway to substantially improve this bound (KATRIN Collaboration, 



J. Angrik et al.|2004 Qtten &: Weinheimer|2008 ) . Cosmological analyses variously 



limit the sum over mass eigenstates to Y2i m v(,i) 0-2 — 0.6 eV ( Abazajian et al 



2011) 



Two types of neutrino mass terms can be added to the SM. Neutrinos can have 
Dirac masses, analogous to those of other SM fermions, if the SM is enlarged to 
include a right-handed neutrino field. Because neutrinos lack charges or other 
additively conserved quantum numbers, lepton-number-violating Majorana mass 
terms can also be added, v^txilvl and v R m R v R , where the former is the only 
dimension-five operator that can be constructed in the SM. (The subscripts L 
and R denote left- and right-hand projections of the neutrino field v, and the 
superscript c denotes charge conjugation.) 



In the seesaw mechanism (Gell-Mann, Ramond Sz Slansky 1979 Mohapatra 



& Senjanovic 1980; Yanagida 1980) the Dirac and Majorana mass terms are 



combined in a manner that provides an attractive explanation for light neutrinos, 



( n ^ 



\ mjj m R J 

where mx ~ in part because of double beta decay constraints. When diagonal- 
ized, the matrix yields heavy and light neutrino mass eigenstates, 



m L ~ me- 



mo 

l 



with the latter related to the typical Dirac mass of the SM by the coefficient 
Trio/mR. If we assume the scale of the new physics which m R represents is 
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S> rriDi then a candidate small parameter mo/mii is available to explain why 
neutrinos are so much lighter than other SM fermions. Small neutrino masses are 
thus explained as a consequence of the scale rriR of new physics beyond the SM. 

Neutrinos of definite mass are the stationary states for free propagation, while 
neutrino flavor eigenstates are produced in weak interactions. Simplifying here 
to two flavors, the relationship of the flavor {z/ e , v^} and mass {v±, v{\ eigenstates 
can be described by a single vacuum mixing angle 9 V , 

u e = cos 6 v \ui) + sin 9 v \u2) = — sin 9 V \ v\) + cos 9 V \ 1/2). (25) 

Consequently, an arbitrary initial state \u(t = 0)) = a e (t = 0)|f e ) + a M (i = 0)|f M ) 
of momentum k ~ E, as it propagates downstream, evolves according to 
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(26) 



where an average overall wave function phase has been removed from the neutrino 
mass matrix (represented in the flavor basis) . For the special case of a v e produced 
at time t = 0, the solution of this equation yields 



Pu e (t) = \{veHt))\ 2 = 1 - sin 2 29 v sin 2 



AE 



1 - - sin 2 20 v (27) 



where the downstream oscillation depends on the difference Srn^i = m\ — m\. 
(If this problem is done properly with wave packets, the oscillation persists until 
the two mass components separate spatially, yielding the asymptotic result on 
the right.) The oscillation length Lq = 47r/ic£ , /5m 2 1 c 4 is shorter than the Earth- 
Sun distance for a typical solar neutrino of energy ~ 1 MeV provided Sm^ > 
1.6 x 1CT 11 eV 2 . Thus solar neutrinos are interesting for oscillation studies because 
of their sensitivity to extremely small neutrino mass differences. 
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4.2 Oscillation Basics: Matter and the MSW Mechanism 



Mikheyev & Smirnov (1985, 1986) showed that the density dependence of the 



neutrino effective mass, a phenomenon first discussed by Wolfenstein ( 1978af b), 



could greatly enhance oscillation probabilities. Their original numerical work 
was soon understood analytically as a consequence of level crossing: a neutrino 
produced in the core as a v e is adiabatically transformed into a by traversing 
a critical solar density where the v e and v„ effective masses cross. It became 
clear that the Sun is not only an excellent neutrino source, but also a natural 
regenerator for enhancing the effects of flavor mixing. 



Equation (26) describing vacuum oscillations is altered in matter 
d ( a e\ if 2E\/2Gfp{x) — brn^y cos 29 v 8m^sxa.20y\ I a e 

dx \a^) AE \5m 2 2l sin 29 v - 2Ey/2G F p{x) + 5m 2 21 cos 26 v ) \ dp 

(28) 

where is the weak coupling constant and p(x) the solar electron number den- 
sity. The new contribution to the difference in diagonal elements, AE^/2Gfp{x), 
represents the effective contribution to m 2 that arises from neutrino-electron scat- 
tering. The indices of refraction of electron and muon neutrinos differ because 
the former scatter via charged and neutral currents, while the latter have only 
neutral current interactions. For 9 V < 7r/4 — the "normal hierarchy" where the 
lighter mass eigenstate makes the larger contribution to v e in vacuum - the mat- 



ter and vacuum contributions to the diagonal elements of Eq. ( 28 ) have opposite 
signs. 



We can diagonalize the right-hand side of Eq. (28) to determine the heavy and 



light local mass eigenstates and eigenvalues m#(x) and i71l(x), functions of p(x) 
\vl(x)) = cos9(x)\u e ) —sin6(x)\u IM ) \vh{x)) = s\n8{x)\v e ) +cos6{x)\v ll ) (29) 
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where the local mixing angle 
■ tar \ sin26 v 



cos 26 (x) 



-X{x) 



(30) 



y / X 2 {x) +sin 2 26 v y / X 2 (x)+sm 2 26 v 

depends on X(x) = 2\[2G ' p p(x)E / 'Sm^ — cos 26 v . Unlike the vacuum case, these 
are stationary states for propagation only if p(x) is constant . Otherwise, defining 



\v(x)) = a H (x)\v H {x)) + a L (x)\u L (x)), Eq. (|28j) becomes 

X(x) ia(x) \ I 

'. ( 31 ) 

-ia(x) — \(x) 

The splitting of the local mass eigenstates and the local oscillation length are 

inkcE 
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l="f 
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I 4E \ 









X(x) = dm^JX 2 ^) + sin 2 26 z 



L (x) 



\{x)c A 



(32) 



while eigenstate mixing is governed by the density gradient 

AE 2 \ V2G F £p(x)sm26 v 



a(x) 



5m 2 21 



(33) 



X 2 (x) + sin 2 2^ 

The splitting achieves its minimum value, 25m 2 l su\26 v , at a critical density 
p c = p{x c ) where X(x) — > 0, 

2V2EG F p c = bm\ x cos 26 v . (34) 



The diagonal elements of the original flavor matrix of Eq. (28) cross at p c . 

The crux of the M SW mechanism is the adiabatic crossing of the critical density, 
illustrated in the left panel of Fig. [6j The adiabatic condition is determined by the 
requirement 7(x) 



am 

a(x) 



S> 1, which allows one to treat Eq. (31) as diagonal 



This condition becomes particularly stringent near the crossing point, 



7c = l(x c ) 



sin 2 26 v 5m 2 1 



1 



cos 26 v 2E 

where H c and L c are the solar density scale height and local oscillation length at 



1 dp(x) I 



27rtan26 , „ > 1, 



(35) 



x c . If the H c > L c , Eq. (31) then yields (Bethe 1986) 



P, 



adiab 



— + - cos 26 v cos 26i 



(36) 
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where Q% = 0{xi) is the local mixing angle at the density where the neutrino is 
produced. The adiabatic solution depends on the local mixing angles where the 
neutrino begins solar core) and ends (9 V , in vacuum) its propagation. 

For illustration, consider the case of a small 6 V ~ 0. A solar v e created in 
the high-density solar core is then nearly identical to the heavy-mass eigenstate 
(Oi ~ 7r/2), provided the vacuum mass difference between the eigenstates is not 



too large (see Eq. (30). If the subsequent propagation is adiabatic, the neutrino 
remains on the heavy-mass trajectory, crossing the critical density (0(x c ) = vr/4), 
and finally exiting the Sun. But in vacuum the heavy mass eigenstate is ~ v^: 
a nearly complete flavor change, v e — > u^, has occurred, through an adiabatic 
rotation of the local oscillation angle from 6i ~ tt/2 to Of = 9 V during propagation. 
If the adiabatic condition is not satisfied, e.g., 7 C < 1, an accurate analytic 



solution can still be obtained (Haxton 1986, Parke 1986). As we have seen, the 
nonadiabatic behavior is governed by the density scale height at x c - One can 
replace the actual solar density by an effective one, e.g., a linear density "wedge" 
that has the correct derivative at x c (thereby incorporating the effects of the 
density gradient at the most sensitive point), while also starting and ending at 
the appropriate initial and final densities (thereby also building in the adiabatic 
limit). The resulting generalization of Eq. ( |36[ ) is 

P»e = \ + \ cos 26 v cos 20; (1 - 2P hop ) P hop = e" 7 ^ 2 (37) 

where Phop; the Landau- Zener probability of hopping from the heavy mass tra- 
jectory to the light trajectory on traversing x c , vanishes in the highly adiabatic 



limit, 7 C » 1 (so that Eq. (J37J) reduces to Eq. (36)). When the crossing becomes 
highly nonadiabatic (7,. < 1 ), Phop — > 1: the neutrino exits the Sun on the light 
mass trajectory, which for small mixing angles means it remains ~ v e . 



41 

Thus strong conversion of solar neutrinos is expected when 1) the propagation 
is adiabatic (7 C > 1) and 2) there is a level crossing (there is enough matter at 
the v e production point that u e (x{) ~ vh{xi))- The right panel of Fig. [6] shows 
the white triangle of parameters in the Srri^i/E — sin 2 26 v plane where both 
constraints are satisfied. Within this triangle, strong conversion occurs. One can 
envision superimposing on this triangle the spectrum of solar neutrinos, plotted 
as a function of dm^i/E for some choice of Srn^i and 6 V . Depending on how that 
spectrum is positioned vertically (a function of Sm^i) or horizontally (a function 
of 6 V ) one can alter the resulting spectrum of u e s in several characteristic ways, 
for example, suppressing the low-energy or high-energy neutrinos preferentially, 
or even (in the case of small mixing angles) those of intermediate energy. 

In early fits to the neutrino data three potential MSW solutions were frequently 
discussed, designated by SMA (small mixing angle: dm^i ~ 5.4 x 10 -6 eV, sin 2 26 ' 
0.006), LMA (large mixing angle: bm\ x ~ 1.8 x 10" 5 eV, sin 2 26 ~ 0.76), and 
LOW (low probability, low mass: bm\ x ~ 7.9 x 10" 8 eV, sin 2 26 ~ 0.96). (The 



parameter values are taken from Bahcall, Krastev & Smirnov ( 1998 ) and are rep- 



resentative of fits done at that time.) These solutions are indicated schematically 
by the colored regions in Fig. |6] The solution consistent with the solar neutrino 
data proved, ironically, to be the LMA solution - not the SMA solution where 
matter effects so greatly enhance the oscillations. 

5 GLOBAL ANALYSES AND NEUTRINO PROPERTIES 

Neutrino physics has made great progress in the last 15 years, as reactor and ac- 
celerator neutrino experiments have added new information to that obtained from 
solar and atmospheric neutrino experiments. The three mixing angles of the 3x3 



42 

neutrino mass matrix, the magnitudes of the two mass differences, and from solar 
neutrino experiments the sign of one of these have all been determined. The phe- 
nomena that can be explored with solar neutrinos were illustrated previously for 
the two-flavor case: flavor oscillations, affected by matter, not only alter fluxes, 
but lead to distinctive spectral distortions, and may produce day-night differences 
due to neutrino passage through the Earth. The various experimental collabo- 
rations as well as independent groups have developed global analysis methods 
to analyze solar neutrino experiments, taking into account the constraints other 
recent measurements have imposed. Here we summarize the conclusions of such 
analyses, relying particularly on work done by the Bari and Valencia groups. 

5.1 Vacuum Mixing Angles and Mass 2 Differences 

In the SM case of three light neutrino flavors, the relationship between flavor 
{v e , u^, v T } and mass { ui, v% } eigenstates is described by the PMNS matrix 
( |Maki, Nakagawa fc Sakata||1962[ |Pontecorvo|l967[ ) 
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(38) 



where Cjj = cos Oij and Sij = sin 9ij . This matrix depends on three mixing angles 
$12) $13; and $23) of which the first and last are the dominant angles for solar 
and atmospheric oscillations, respectively; a Dirac phase S that can induce CP- 
violating differences in the oscillation probabilities for conjugate channels such 
as — > v e versus — > u e ; and two Majorana phases a\ and 02 that will affect 



the interference among mass eigenstates in the effective neutrino mass probed in 
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the lepton-number- violating process of neutrinoless double (3 decay. 

It became apparent in early analyses that combined solar and reactor neutrino 
data in two-flavor analyses that there was some hint of the third flavor, a nonzero 
#13. The KamLAND Collaboration analysis employed the three- flavor v e survival 



probability of Fogli et al. (2000) in which the influence of #13 in modifying the 



two- flavor result is explicit, 



Pf e = COS 4 013 Pee + Sm 4 9 13 



2u 



(39) 



where P^ is the two-flavor survival probability in matter evaluated for the 
modified electron density p(x) — > p{x) cos 2 #13. The analysis yielded sin 2 #13 = 



0.020 ± 0.016 (KamLAND Collaboration, A. Gando et al. 2011), a result con- 



sistent with the long-baseline z/ e appearance results announced shortly after- 



wards, 0.008 < sin 2 #13 < 0.094 (T2K Collaboration, K. Abe et al. 2011) and 



0.003 < sin 2 0i3 < 0.038 (MINOS Collaboration, P. Adamson et al. 2011). In 



2012 results from reactor v e disappearance experiments became available, yield- 
ing sin 2 0i3 = 0.022 ± O.Oll(stat) ± 0.008(sys) (|Double Chooz Collaboration, Y. 



Abe et al.|2012| , sin 2 0i 3 = 0.0236 ± 0.0042(stat) ±0.0013(sys) ( |Daya Bay Collab- 
oration, F. P. An et al.||2012| ), and sin 2 13 = 0.0291 ± 0.0035(stat) ± 0.0051 (sys) 



(RENO Collaboration, J. K. Ahn et al. 2012). The latter two results, because 



of their precision, effectively remove a degree of freedom from three-flavor solar 
neutrino analyses. 

The mass differences and mixing angles from the global analyses of the Bari 
( Fogli et al.||2012[ ) and Valencia ( |Forero, Tortola fc VaJle|2012| | groups, including 
experimental results through the Neutrino 2012 Conference, are shown in Table 
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[3j The two analyses are generally in quite good agreement and yield (in degrees) 

(40) 



'12 ~ \ 



33.61^ 



713 ~ \ 



34.41J;? 



8.96^? 



9.06j£|? 



(41) 



Bari 
Valencia 

The agreement in the solar neutrino mass difference Sm^i is also excellent, 

(7.54±g;||) x 10- 5 eV 2 Bari 
(7.62ig;Jg) x 10~ 5 eV 2 Valencia 
The values for 9±2 and Sm^ he in the LMA region of Fig. [6j Srri^i corresponds, 
for 10 MeV neutrinos, to an MSW crossing density of ~ 20 g/cm 3 , or equivalently 
a solar radius of r ~ O.24i? , the outer edge of the Sun's energy-producing core. 
The crossing density for the atmospheric Sm^, again for 10 MeV neutrinos, is 
~ 1.6 x 10 3 g/cm 3 . Thus this crossing requires electron densities far beyond those 
available in the Sun - though typical of the carbon zone in the mantle of a Type 
II supernova, where this second crossing plays a significant role. 

These global analysis results can be compared with those from the recent SNO 
three-flavor combined analysis, which used all available solar neutrino data and 
the results from KamLAND. This analysis, summarized in Fig. [7| gives at la 

sin 2 #12 = 0.308 ± 0.014 5m 2 21 = (7.41+°;^) x 10" 5 eV 2 sin 2 13 = 0.025±g;g^| 

(42) 

These values are in excellent agreement with the corresponding la Bari and 
Valencia results of Table [3j the SNO combined analysis and Bari best values 
match particularly well. The main consequence of the inclusion of new reactor 
and accelerator results in the global analyses is a substantial reduction in the 
uncertainty on #13. 
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5.2 Spectral Distortions: LET Analyses and Borexino 

Characteristic spectral distortions are one of the signatures of oscillations in mat- 



ter. Rather fortuitously, if one evaluates Eq. (34) for the neutrino energy where 



the MSW critical density corresponds to the electron density at the center of the 
Sun, p ~ 6 x 10 25 /cm 3 , one finds ££rit ~ ^ MeV, an energy in the center of the 
solar neutrino spectrum. Neutrinos below this energy will not experience a level 
crossing on exiting the Sun, and thus will oscillate approximately as they would 
in vacuum - an average (two-flavor) survival probability of 

pvacuum ~ l _ I s i n 2 26 12 ~ 0.57, (43) 

using #i2 ~ 34°. This can be compared to the matter-dominated survival proba- 
bility, appropriate for neutrinos much above the critical energy, 

phigh density ^ ^2 ^ _ Q ^ ^ 

Most of the 8 B neutrinos studied by SNO and Super-Kamiokande undergo matter- 
enhanced oscillations. The matter/vacuum transition predicted by the MSW 
mechanism can be verified by comparing the survival probabilities of low-energy 
(pp or 7 Be) and high-energy ( 8 B) neutrinos. Alternatively, if the thresholds in 
SNO and Super-Kamiokande are lowered sufficiently, spectral distortions will be 
detectable in the 8 B spectrum: low-energy 8 B neutrinos coming from the outer 
core will not experience a crossing, and thus will have a higher survival rate. 

The flux of low-energy pp neutrinos is well constrained in global analyses be- 
cause these neutrinos dominate the SAGE and GALLEX/GNO counting rates. 
(The need for an elevated survival probability for these neutrinos was an im- 
portant factor in early model-independent analyses that concluded undistorted 
neutrino fluxes could not account for the data.) Furthermore Borexino has now 
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provided a direct, exclusive measurement at a precise energy, corresponding to 
the 860 keV neutrinos from 7 Be electron capture. 

To probe lower energy 8 B neutrinos, a joint re-analysis of Phase I and Phase 
II data from the Sudbury Neutrino Observatory was carried out with an effective 



kinetic energy threshold of T c r = 3.5 MeV (SNO Collaboration, B. Aharmim 



et al. 2010). While the low-energy threshold analysis (LETA) had several mo- 
tivations (e.g., the enlarged data set improved the overall precision of the flux 
determinations), a principal goal was enhancing prospects for detection of the 
predicted upturn in P Ve with decreasing neutrino energy. An effort similar to the 
LETA analysis is now underway in Super-Kamiokande IV: preliminary results 



were recently described by Smy (2012) 



Figure [8] summarizes the data. The pattern defined by the pp v e flux deduced 
from global analyses, the 7 Be v e flux derived from the Borexino ES measurement, 
and the SNO results are generally in good agreement with the expected MSW 
survival probability. However, while the SNO LETA band is compatible with the 
MSW prediction, the band's centroid trends away from the theory with decreasing 
energy. 

5.3 Day-Night Differences 

Two sources of time variation in neutrino rates are the annual ~ 7% modulation 
associated with the 1.7% eccentricity in the Earth's orbit around the Sun, and 
the daily variation associated with terrestrial matter effects, which influence the 
night-time flux of up-going neutrinos. Both effects have been the subject of careful 



experimental studies: see, e.g., SNO Collaboration, B. Aharmim et al. (2005); 



Super-Kamiokande Collaboration, M. B. Smy et al. (2004). The integrated day 
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night asymmetry in neutrino detection rates 

where i?o and Rj$ denote the day and night rates, is the quantity most often 
studied to assess matter effects associated with solar neutrino passage through the 
Earth. In principle, similar differential quantities could be defined as functions 
of the neutrino energy and zenith angle. However, the detection of even the 
integrated difference ^4dn is statistically challenging, as the effect is expected to 
be only a few percent. 

^4dn provides an "on-off" test where the matter effects can be measured di- 
rectly, unlike the solar case where matter effects must be deduced from phenom- 
ena such as spectral distortions. The magnitude of the neutrino regeneration 
associated with passage through the Earth depends on the neutrino energy, the 
assumed oscillation parameters, and, to some extent, detector location, as that 
determines the possible trajectories through the Earth to the Sun. 

The high counting rate of Super-Kamiokande is an advantage in constraining 
the contribution to Adn from the high-energy 8 B neutrinos. The expected Super- 
Kamiokande ES ^4dn is illustrated in Fig. [9] as a function of Snv^x and sin 2 #13, 
in the two- and three-flavor cases. For the current global best-fit parameters, the 
effect should be about —3%. The results for the four phases 

-2.1 ±2.0 ±1.3% SKI 6.3 ±4.2 ±3.7% SK II 

-5.9 ±3.4 ±1.3% SKIII -5.2 ±2.3 ±1.4% SK IV (preliminary) 

(46) 

yield a combined result of Agf,(SK) = -4.0 ± 1.3 ± 0.8, in good agreement with 



expectations but still consistent with no effect at 2.6<7 (Smy 2012). 



The SNO Collaboration has analyzed day-night effects in the v e channel in 



48 

their combined analysis, approximating the effect as linear in E v , 

^ N (SNO) = -a -a 1 ( I ^ v -10). (47) 

The null hypothesis that there are no day/night effects influencing the v e survival 
probability (so ao = 0, a\ = 0) yielded a Ax 2 = 1.87 (61% c.l.) compared to the 



best fit (SNO Collaboration, B. Aharmim et al. 2012). The Borexino result for 



the ES at 862 keV flBorexino Collaboration, G. Bellini et al.||2012b| ) , 



^dn(862 keV) = -(0.001 ± 0.012 ± 0.007), (48) 

is consistent with the expectation that | Ag|j(862 keV| < 0.001 for the LMA 
solution at this energy. 

6 NEUTRINO CONSTRAINTS ON SOLAR STRUCTURE 

One of the important consequences of the increasingly precise understanding of 
neutrino flavor physics is the opportunity to return to one of the early goals of 
solar neutrino spectroscopy, using the neutrino as a probe of the physics of the 
solar interior. Neutrino fluxes, sensitive to nuclear reaction rates and core tem- 
perature, can be combined with helioseismic observations, sensitive to radiative 
opacities and microscopic diffusion, to place stringent constraints on the SSM 
and to test some of its implicit assumptions. This program is of broad signifi- 
cance to stellar astrophysics because the SSM is a particular application of the 
general theory of main-sequence stellar evolution. Because we know the Sun's 
properties far better than those of any other star, the SSM provides one of our 
best opportunities to test that theory against precise data, and thus to identify 
shortcomings. 

A decade ago the SSM was in spectacular agreement with observations apart 
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from solar neutrino data, a fact that supported suggestions that the solar neutrino 
problem might have a nonsolar origin. However, as the SSM makes a number 
of simplifying assumptions, it is perhaps inevitable that some experimental test 
of our Sun will eventually demonstrate the model's shortcomings. Over the past 
decade, the development of 3D hydrodynamic models of near-surface solar con- 
vection, a more careful selection of spectral lines, and, in some cases, relaxation 
of the assumption of local thermodynamic equilibrium in line formation have led 
to significant changes in the analysis of data on photospheric absorption lines. 
The most recent revisions reduced the abundances of the volatile CNO elements 
and Ne by ~ 0.10-0.15 dex (equivalently, by ~ 25 — 40%), relative to older compi- 
lations of solar abundances - though significant debate continues. The differences 
between the new AGSS09 and the older GS98 abundances can be summarized in 
the respective total-metal-to-hydrogen ratios of (Z/X)© ~ 0.018 and 0.023, re- 
spectively. As the SSM assumes a homogeneous zero-age composition, adoption of 
the AGSS09 abundances produces a modern Sun with a lower core metallicity, af- 
fecting solar neutrino flux predictions and substantially degrading the agreement 
between the SSM sound velocity profile and that deduced from helioseismology. 

6.1 The Solar Abundance Problem and its SSM Implications 

Past studies of SSMs with low Z/X interiors similar to that of the AGSS09- 
SFII SSM have revealed a number of difficulties with respect to observation: 
the radius of the convective zone boundary Rcz, the interior sound speed and 
density profiles, and the surface He abundance Y5 all move outside the ranges 



determined from helioseismic analyses (Bahcall, Serenelli & Pinsonneault 2004 



Basu Antia|20"04 Montalban et al.|2004 ). These changes reflect the increase in 
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the efficiency of radiative transport and decrease in core molecular weight found 
in low-Z models. For example, Rcz moves outward in low-Z models because 
radiative transport dominates over a larger fraction of the solar interior. Similarly, 
Ys decreases: As SSM energy generation is fixed by the measured luminosity 
Lq, the cooler core in low-Z models must be compensated by an increase in the 
available fuel X, and consequently a lower core Y and thus surface Ys- In contrast, 
SSM predictions using the older, high-Z GS98 abundances are in much better 
agreement with observation. (See Tableland Fig. [5|) The inconsistency between 
the SSM parameterized using the best current description of the photosphere 
(AGSS09 abundances) and the SSM parameterized to optimize agreement with 
helioseismic data sensitive to interior composition (GS98 abundances) is known 
as the solar abundance problem. 

The solar abundance problem could have a pedestrian solution: the 3D analysis 



of |Caffau et al.| ((2008} |2009[ |2010| ) yielded abundances higher than the AGSS09 
values, though this appears to be due to spectral line choices rather than photo- 
spheric model differences ( Grevesse et al.pOlT ) , a conclusion supported by a re- 



cent comparison between solar model atmospheres computed by different groups 



(Beeck et al. 2012). Alternatively, an upward adjustment in associated atomic 
opacities could compensate for a low-Z interior, if some justification for such a 
change could be identified. 

However, the solar abundance problem could be more fundamental. An im- 
portant assumption of the SSM - a homogeneous zero-age Sun - is not based on 
observation, but instead on the theoretical argument that the proto-Sun likely 
passed through a fully convective Hayashi phase as the pre-solar gas cloud col- 
lapsed, thereby destroying any composition inhomogeneities that might have ex- 
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isted. Yet we know that chemical inhomogeneities were re-established during 
solar system formation: processes operating in the proto-planetary disk removed 
~ 40 — 90M® of metal from the gas, incorporating this material in the gaseous gi- 
ants ( Guillot|2005 ). The gas from which these metals were scoured - representing 



perhaps the last ~ 5% of that remaining in the disk - would have been depleted 
in metal, and enriched in H and He. The fate of that gas is unknown, but if 
it were accreted onto the Sun, it plausibly could have altered the composition 
of the convective zone, depending on the timing of the accretion and thus the 
maturity of the proto-Sun's growing, chemically segregated radiative core. That 
is, as there is a candidate mechanism for altering the convective zone late in 
proto-solar evolution, involving enough metal to account for the AGSS09/GS98 
differences, it is not obvious that the SSM assumption of homogeneity is correct. 
The solar abundance problem has three connections to neutrinos: 

1. Neutrino fluxes are sensitive to metallicity, and thus can be used to cross- 
check the conclusions drawn from helioseismology. Below we describe what 
the neutrino fluxes currently tell us. 

2. Neutrino fluxes place important constraints on nonstandard solar models 
(NSSMs) motivated by the solar abundance problem, such as those recently 
developed to explore accretion from the proto-planetary disk. 

3. Planned measurements ( Chen|[2~006 Franco|[2011 ) of the CN solar neutrino 



flux have the potential to directly measure the solar core abundance of C+N 
with a precision that will impact the solar abundance problem. 

In principle the temperature-dependent 8 B and 7 Be neutrino fluxes have suffi- 
cient sensitivity to metallicity to impact the solar abundance debate. As Table [2] 
shows, the AGSS09-SFII and GS98-SFII SSMs differ by 21.6% and 9.6% in their 
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8 B and 7 Be flux predictions, respectively, which one can compare to the 14% and 
7% SSM uncertainties on these fluxes obtained by varying SSM input parameters 
according to their assigned errors. These total SSM uncertainties were deter- 
mined by adding in quadrature the individual uncertainties from 19 SSM input 
parameters, including abundance uncertainties as given in the respective solar 
abundance compilations. In the case of the 8 B flux, the important uncertainties 
include those for the atomic opacities (6.9%), the diffusion coefficient (4%), the 
nuclear S-factors for 3 He+ 4 He (5.4%) and 7 Be+p (7.5%), and the Fe abundance 
(5.8%). 

Unfortunately the current neutrino data do not favor either abundance set. 
The last entries in Table [2] give the x 2 functions and compatibility functions 



_P agr obtained in the SSM for the two sets, following Gonzales-Garcia, Maltoni 



& Salvado (2010) but using the updated Solar Fusion II S-factors. The two 
models are identical in quality of fit to the data, with XAGSS09-SFII = 3-4 and 
^GS98-SFii = 3-5. From Table [2] one sees that a SSM model intermediate in 
metallicity between AGSS09 and GS98 would optimize the agreement. 

6.2 Solar Models with Accretion 



Guzik (2006), Castro, Vauclair & Richard (2007), and Guzik k Mussack (2010) 



considered the possibility that the solar abundance problem might be due to 
accretion of metal-poor matter onto the Sun's convective envelope, diluting its 



pre-solar composition. Haxton & Serenelli (2008) suggested the mechanism for 



such dilution described above, accretion of disk gas from which the metals were 
previously scoured in the process of planetary formation. Evidence supporting 
the hypothesis that planet formation could affect the surface metallicity of the 
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host star was offered by Mefendez et al. (2009), who found that the peculiar 



differences in the surface abundances of the Sun, measured relative to similar 
stars (solar twins), correlate with condensation temperatures, and thus plausibly 



with disk chemistry. Nordlund (2009) argued this accretion scenario might also 
provide a natural explanation for the anomalous metallicity of Jupiter. 

While the process of planetary formation is not well understood, the standard 
picture invokes a chemically differentiated thin disk, with dust, ice, and thus 
metals concentrated in the midplane, and outer surfaces dominated by H and 
He. This configuration arises from a combination of gas cloud collapse that is 
rapid at the poles but inhibited by angular momentum at the equator, and gas 
cooling that allows differential condensation of various elements as ice or grains 
according to their condensation temperatures. The formation of planetesimals in 
the midplane and their self-interactions lead to the formation of rocky cores of 
planets. The gas giants, which are sufficiently distant from the Sun that ice can 
augment their cores, reach masses where tidal accretion of gas can further feed 
their growth. 

The implications of this disk chemistry for solar initial conditions is difficult to 
assess because of a number of uncertain parameters. While we know the planets 
are substantial metal reservoirs, we do not know the fate of the depleted gas that 
dominates the disk surface: it might have been removed by the solar wind, or 
alternatively deposited on the Sun through magnetospheric mechanisms, such as 
those operating in young accreting T Tauri stars. If accretion occurs, its timing 
relative to early solar evolution is critical. The timescale for planetary formation 
is generally estimated at 1-10 M.y, while the SSM predicts that the early Sun's 
convective boundary moves outward, in response to the growing radiative core, 
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over a longer period of tqz ~ 30 M.y. If accretion occurs early in this period, when 
most of the Sun's mass is in its convective envelope, any resulting nonuniformity 
in solar abundances would be negligible. But if accretion occurs later when 
the convective zone is similar to that of the modern Sun, and thus contains 
only ~ 2% of the Sun's mass, the chemical processing of ~ 0.05M Q of gas, 
dust, and ice in the planetary nebula could have very significant consequences. 
Alternatively, as suggested by hydrodynamic simulations of the contraction of 
the protosolar nebula ( Wuchterl Klessen|2001 ) and stellar models that include 
episodic accretion during the pre-main-sequence phase with large mass accretion 



rates (Baraffe & Chabrier 2010), the Sun may have avoided the fully convective 
phase altogether. Finally, the mass and composition of the accreted material 
is highly uncertain. Because condensation temperatures for the elements vary 
widely - e.g., from ~1400 K for Fe to ~300 K for C - the composition would 
likely evolve with time as the disk cools. 

Motivated in part by such considerations, a NSSM was recently developed to 
test whether early accretion from a planetary disk could resolve the solar abun- 



dance problem (Serenelli, Haxton & Peha-Garay 2011). The work illustrates the 
importance of neutrino physics in limiting NSSMs. In SSMs the pre-solar compo- 
sition parameters Y; n i and Z m \ as well as the mixing length parameter omlt are 
determined by iterating the model to reproduce the present-day solar luminos- 
ity Lq, radius Rq, and surface metal-to- hydrogen ratio Z5/X5. This algorithm 
was modified for the accreting NSSM by allowing for a mass M ac of accreted 
material with fixed composition (X ac , Y ac , Z ac ), deposited on the early Sun uni- 
formly, beginning at time r aCi i and lasting a time Ar ac . Prior to time r aCj i the Sun 
was evolved as a SSM with mass M — M ac and composition defined by (Xj n i, 
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Ymi, Zi n i). During the subsequent accretion phase, the simplifying assumption 
Xac/Yac = Xi n i/Yi n i was made. For a given set of fixed accretion parameters, 
Ymi, Zj n i, and «mlt were then adjusted iteratively in order to reproduce L , 
Rq, and Zs/Xs, as in the SSM. Thus Y^, Zj n i, and «mlt become functions of 
the assumed accretion parameters, so that in place of a single SSM solution, a 
family of solutions is obtained with different interior compositions. Unphysical 
solutions are rejected, e.g., the AGSS09 Z5/X5 is not compatible with the accre- 
tion of large quantities of metal-rich material at late times, when the convective 
envelope is thin. 



Fig. 10 shows the helioseismic consequences of late accretion - accretion onto 
a young Sun with a thin convective envelope similar to that of the modern Sun. 
Also shown are the variations in Yi n ; and Zi n j, the initial core He and heavy- 
element mass fractions, that can be achieved through accretion. The helioseismic 
observables are very constraining. The candidate accretion solution to the solar 
abundance problem that one might naively envision - a low-Z surface consistent 
with AGSS09 and an interior similar to GS98, with higher Zj n i ~ 0.019 and Yi n j ~ 
0.28 - can be achieved with metal-free and metal-poor accretion involving modest 
accreted masses M ac ~ O.OIMq. These models bring Y$ into good agreement 
with observation and produce some improvement in the sound-speed figure-of- 
merit {5c/ c). But the lower Zg that accompanies metal- free/metal-poor accretion 
allows the convective zone radius to move outward, exacerbating the existing 
AGSS09-SFII SSM helioseismic discrepancy in R$. 

Neutrino flux measurements impose a second class of constraints on accreting 
NSSMs. The scenario discussed above - metal-free or metal-poor accretion with 
Mac ~ O.OIM0 - can marginally improve the agreement with experiment, con- 
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sistent with earlier observations that this produces an interior similar to GS98, 
and that the GS98 and AGSS09 SSMs yield equivalent fits to the solar neutrino 
data. But with larger M ac the agreement quickly deteriorates as the resulting 
high-Z interior leads to rapid increases in the 8 B and 7 Be neutrino fluxes. The 
power of contemporary neutrino flux measurements to constrain NSSMs is quite 
remarkable. Large classes of accretion parameters - mass, time, composition, 
and duration - lead to modern Suns with the proper luminosity and radius and 
the AGS009 Zg, while satisfying the underlying equations of stellar structure. 



Yet very few of these solutions produce acceptable neutrino fluxes, as Fig. 11 
illustrates. 

Effectively the recent progress made on neutrino mixing angles and mass dif- 
ferences has made the neutrino into a well understood probe of the Sun. We now 
have two precise tools, helioseismology and neutrinos, that can be used to see into 
the solar interior, complementing the more traditional astronomy of solar surface 
observations. Effectively we have come full circle: the Homestake experiment 
was to have been a measurement of the solar core temperature, until the solar 
neutrino problem intervened. 

6.3 Neutrinos as a Direct Probe of Solar Composition: SNO+ 

While solar fusion is dominated by the pp chain, the CN cycle generates ~ 1% 
of solar energy as well as 13 N and 15 neutrino fluxes of 2.96 (2.17) and 2.23 
(1.56) xl0 8 /cm 2 s, respectively, for the GS98-SFII (AGSS09-SFII) SSM. These 
fluxes can be measured in scintillation detectors in an energy window of 800-1400 
keV, provided backgrounds in the detectors are reduced to acceptable levels, in- 
cluding in particular interference from the decay of 11 C produced by penetrating 
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muons. Borexino has already established a limit on the CN fluxes, and is pur- 
suing strategies to improve the measurement by vetoing interfering backgrounds 



(Franco 2011 ). A new scintillation detector under construction in the cavity previ- 



ously occupied by SNO will have three times the volume of Borexino and 1 /70th 



the muon background, due to SNOLab's two-kilometer depth (Chen 2006). A 
SNO+ CN neutrino signal/ 11 C background of ~ 10 is expected. 

A measurement of the CN neutrinos would test our understanding of hydrogen 
fusion as it occurs in main-sequence stars substantially more massive than the 
Sun. It could also play an important role in the solar abundance problem. 

The CN neutrino fluxes, like the 8 B and 7 Be neutrino fluxes, depend sensitively 
on core temperature Tq and thus respond to variations in SSM input parame- 
ters that affect Tq- But in addition, the CN fluxes have a linear dependence 
on the pre-solar core abundances of C and N that is unrelated to Tq, reflecting 
the proportionality of CN neutrino fluxes to the abundances that catalyze the 
hydrogen burning. The Tc power-law relationships for the CN and 8 B neutrino 
fluxes and the additional linear dependence of the CN neutrino flux on the abun- 
dance of C+N can be exploited to relate solar neutrino measurements to the 



Sun's pre-solar C and N abundances (Haxton Sz Serenelli 2008) 



ri5 



O) 



15Q)SSM 



»B) 



0(»B) 



SSM 



0.729 



XC+N 



x [1 ± 0.006(solar) ± 0.027(D) ± 0.099(nucl) ± 0.032(6»i 2 )] (49) 

where xq+n is the C+N number abundance normalized to its SSM value. The 
uncertainties were derived from SSM logarithmic derivatives, as described in Sec. 
[2] The first two of these represent variations in all SSM parameters other than the 
nuclear cross sections - including Lq, the opacity, solar age, and all abundances 
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other than C and N, using abundance uncertainty intervals of 
. | Abundance, 0598 - Abundancef GSS09 



(50) 



(Abundancef S98 + Abundancef GSS09 )/2 
Apart from the diffusion parameter D, the net effect of the variations in these 
quantities is an uncertainty of 0.6%: we have formed a ratio of fluxes that is 
effectively insensitive to Tq. The diffusion parameter D is an exception because 
our expression relates contemporary neutrino flux measurements to the pre-solar 
number densities of C and N, and thus must be corrected for the effects of diffusion 
over 4.6 Gy. The differential effects of diffusion on the ratio creates an uncertainty 
of 2.7%, the only significant nonnuclear solar uncertainty. 



Equation (49) is written for instantaneous fluxes, and thus must be corrected 
for the energy-dependent effects of oscillations. The SNO combined analysis 
result 0i2 = Umtlit - or equivalently the Bari or Valencia global analysis 
results of Sec. [5]- imply a < 3.2% uncertainty in the flux comparison of Eq. (49). 
Finally, there are nuclear physics uncertainties. These dominate the overall error 
budget, with the combined (in quadrature) error reflecting a 7.2% uncertainty 
from the 1 N(p,7) reaction and a 5.5% uncertainty from 7 Be(p,7). 

The SNO and Super-Kamiokande measurements of the 8 B flux have reached 
a precision of 3%. SNO+ has the potential to measure the 15 flux to about 
7% in three years of running. Assuming such a result from SNO+ and combin- 
ing all errors in quadrature, one finds that the pre-solar C+N abundance can 
be determined to ± 13%. The precision could be improved substantially by ad- 
dressing the uncertainties in the S-factors S114 and S17. But even without such 
improvements, the envisioned SNO+ result would have a major impact, given 
the existing ~ 30% differences in the AGSS09 and GS98 C and N abundances. 
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7 OUTLOOK AND CHALLENGES 

In this review we have summarized the results of 50 years of work on solar neu- 
trinos. The field has been characterized by very difficult experiments carried out 
with great success, producing results fundamental to two standard models, our 
standard theory of stellar evolution and our standard model of particle physics. 
Thirty years of debate over the origin of the solar neutrino problem - a misunder- 
standing of the structure of the Sun, or an incomplete description of the properties 
of the neutrino - ended with the discovery of massive neutrinos, flavor mixing, 
and MSW distortions of the solar neutrino spectrum. The quest to resolve the 
solar neutrino problem spurred the development of a new generation of active 
detectors of unprecedented volume and radiopurity - SNO, Super-Kamiokande, 
and Borexino - that have made solar neutrino spectroscopy a precise science. 
This technology has opened up new possibilities. 

The program of solar neutrino studies envisioned by Davis and Bahcall has 
been only partially completed. Borexino has extended precision measurements 
to low-energy solar neutrinos, determining the flux of Be neutrinos to 5%, and 
thereby confirming the expected increase in the v e survival probability for neu- 
trino energies in the vacuum-dominated region. First results on the pep neutrino 
flux have been obtained, as well as a limit on the CN neutrino fluxes. But a 
larger, deeper version of Borexino, SNO+, will likely be needed to map out the 
low-energy solar neutrino spectrum in detail, including the CN neutrino contribu- 
tions. The observation of these neutrinos in the Sun would provide an important 
test of the nuclear reactions we believe dominate energy generation in massive 
hydrogen-burning stars. 

There are challenging tasks remaining. The flux of solar pp neutrinos, known 
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theoretically to ± 1%, is our most accurately predicted source of v e s. The detec- 
tion of these neutrinos ( Grieb et al.|2011 ) would check the luminosity condition so 
widely used in solar neutrino analyses, the equivalence of solar energy production 
as measured in emitted photons to that deduced from the neutrinos - a test also 
relevant to possible new physics, such as sterile neutrinos (Abazajian et al.|2012p . 
The pp neutrinos could also provide our most precise low-energy normalization 
for calibrating the MSW effects that turn on with increasing energy. We have 
discussed the excellent prospects that SNO+, by measuring the 15 CN neutrinos 
(1.73 MeV endpoint), will directly determine the pre-solar C+N content of the 
solar core to a precision that will impact the solar abundance problem. Though a 
much more challenging task, the measurement of the lower energy 13 N neutrinos 
(1.20 MeV endpoint) to high precision could yield separate determinations of the 
core C and N abundances: the ~ 35% higher flux of 13 N neutrinos reflects the 
ongoing burning of pre-solar carbon in the Sun's cooler outer core, as discussed 
in Sec. [2] Finally, there is the next-to- next-generation challenge of measuring the 



Sun's thermal neutrinos (Haxton & Lin 2000), a solar source that becomes sig- 
nificant below 10 keV, with a peak flux density of ~ 10 9 /cm 2 /s/MeV. Neutrinos 
of all flavors are produced by various Zq — > vv processes operating in the Sun. 
Their fluxes are a sensitive thermometer for the solar core, independent of nuclear 
physics, and depend on the core abundance of Fe and other heavy elements, due 
to the contribution of free-bound transitions where a Zq is radiated. The most 
likely candidate for detection may be the z^ e s, as the correspondence between 
thermal neutrino energies and typical inverse atomic scales suggests some form 
of resonant absorption process. 

In retrospect, it is remarkable that we have learned so much fundamental 
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physics from an object as manifestly complicated as the Sun - with all of its 
3D physics, magnetic fields, convective and radiative transport, and nontrivial 
nuclear and atomic physics. Apart from improving the precision of input pa- 
rameters and small tweaks in the physics, such as the inclusion of He and heavy 
element diffusion, the SSM that helped us extract new neutrino physics from the 
Sun is the same model Bahcall, Fowler, Iben, and Sears employed in 1963 to 
make the first prediction of solar neutrino fluxes. It would not be surprising to 
see cracks finally appearing in that model, under the relentless pressure of ex- 
periments measuring neutrino fluxes with few percent precision and helioseismic 
mappings of interior sound speeds to a few parts in a thousand. 

We have reviewed one possible crack in some detail - the solar abundance 
problem. Despite the lack of consensus on this issue, it has raised an interesting 
question: why have we done so well with a model that assumes a homogeneous 
zero-age Sun, when the large-scale segregation of metals is apparent from plane- 
tary structure? The abundance problem has been linked to dynamical issues in 
the collapse of the solar nebula, the metallicity of solar twins, and the anomalous 
composition of Jupiter. We also have seen that neutrinos could provide a crucial 
check on the possibility of an inhomogeneous Sun because of the sensitivity of 
the CN neutrino flux to the C+N content of the core. It could well turn out that 
we are encountering the first evidence that a more complete model is needed - 
perhaps one that moves beyond an isolated Sun, to address the chemical and dy- 
namical coupling of the Sun and planets at their births. Just as the SSM allowed 
us to exploit our best understood star to test the general theory of main-sequence 
stellar evolution, an effort to develop a Standard Solar System Model - one with 
gas cloud collapse, disk formation and accretion, realistic pre-solar evolution, 
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the growth of planets, and the coupled astrochemistry of the Sun and planets - 
could provide a needed template for interpreting what we are now learning about 
exoplanets and their host stars. 
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Table 1: SSM characteristics are compared to helioseismic values, as determined 



by Basu & Antia (1997, 2004). X, Y, and Z are the mass fractions in H, He, and 
metals. The subscripts S, C, and ini denote current photospheric, current core, 
and zero-age values. Rcz is the radius to the convective zone, and {Sc/c) is the 
average fractional discrepancy in the sound speed, relative to helioseismic values. 



Property 


GS98 


AGSS09 


Solar 


(z/x) s 


0.0229 


0.0178 




Z.s- 


0.0170 


0.0134 




Ys 


0.2429 


0.2319 


0.2485±0.0035 


Rcz/Rq 


0.7124 


0.7231 


0.713±0.001 


(8c/c) 


0.0009 


0.0037 


0.0 


1c 


0.0200 


0.0159 




Yc 


0.6333 


0.6222 




Zini 


0.0187 


0.0149 




Y- ■ 

1 mi 


0.2724 


0.2620 
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Figure 1: (Color online) The left frame shows the three principal cycles com- 
prising the pp chain (ppi, ppll, and ppIII), the associated neutrinos that "tag" 
each of the three branches, and the theoretical branching percentages defining 
the relative rates of competing reactions (GS98-SFII SSM). Also shown is the 
minor branch 3 He+p — > 4 He+e + +^ e , which generates the most energetic neu- 
trinos. The right frame shows the CN I cycle, which produces the 13 N and 15 
neutrinos. 
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Figure 2: (Color online) Left panel: The data, the best quadratic+screening 
result for 833(E), and the deduced best quadratic fit (red line) and allowed range 
(yellow band) for Sl^ 1 " 6 . Right panel: Sn(E) vs. center-of-mass energy E, for 
E < 1250 keV. Data points are shown with total errors, including systematic 
errors. The dashed line is based on a theoretical calculation, scaled to fit the 
data. See Solar Fusion II for references and details. 



Table 2: SSM neutrino fluxes from the GS98-SFII and AGSS09-SFII SSMs, with 
associated uncertainties (averaging over asymmetric uncertainties). The solar 
values come from a luminosity-constrained analysis of all available data by the 
Borexino Collaboration. 



v flux 


E* lax (MeV) 


GS98-SFII 


AGSS09-SFII 


Solar 


units 


p+p-> 2 H+e++v 


0.42 


5.98(1 ±0.006) 


6.03(1 ±0.006) 


6-05(1™) 


10 10 /cm 2 s 


p+e~+p— » 2 H+v 
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1.44(1 ±0.012) 


1.47(1 ±0.012) 


i-46(li°;P 


10 8 /cm 2 s 


7 Be+e~-> 7 Li+y 


0.86 (90%) 
0.38 (10%) 


5.00(1 ±0.07) 


4.56(1 ±0.07) 




10 9 /cm 2 s 


8 B-> 8 Be+e++^ 


~ 15 


5.58(1 ±0.14) 


4.59(1 ±0.14) 


5.00(1 ±0.03) 


10 6 /cm 2 s 


3 He±p-> 4 He±e+ +u 


18.77 


8.04(1 ±0.30) 


8.31(1 ±0.30) 




10 3 /cm 2 s 


13 N-> 13 C±e++^ 


1.20 


2.96(1 ±0.14) 


2.17(1 ±0.14) 


< 6.7 


10 8 /cm 2 s 


15 0-> 15 N+e++^ 


1.73 


2.23(1 ±0.15) 


1.56(1 ±0.15) 


< 3.2 


10 8 /cm 2 s 


17 F^ 17 0±e++y 


1.74 


5.52(1 ±0.17) 


3.40(1 ±0.16) 


< 59. 


10 6 /cm 2 s 


X 2 /P agr 




3.5/90% 
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Table 3: Results from global 3f analyses including data through Neutrino2012. 





Bari 
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Fogli et al. 


2012 
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Forero, Tortola & Valle 


2012 




Parameter/hierarchy 


Best la Fit 
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Best 1(7 Fit 
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7 54+0.26 


7.15 o 8.00 


6.99 o 8.18 
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7.27 o 8.01 


7.12 o 8.20 




Sm 2 31 (W- 3 eV 2 


) NH 
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o cr:+0.06 
Z - OO -0.09 
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Figure 3: (Color online) The solar neutrino spectrum, along with the SSM un- 



certainties (Serenelli, Haxton Sz Peha-Garay 2011). A weak branch from the /? 



decay of 1 F that contributes from the CN II cycle is included. The units for the 
continuous sources are cm -2 s _1 MeV _1 . 
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Figure 4: (Color online) Comparison of the measured neutrino rates for the 
chlorine, Kamioka II /III, and SAGE/GALLEX/GNO experiments with the con- 



temporary SSM GS98-SFII, assuming unoscillated fluxes (Serenelli, Haxton & 



Peha-Garay 2011). 
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Figure 5: (Color online) The relative sound speed (c(r) so i ar — c(r)ssivi)/c(r)ssM 
where c(r)ssM is the SSM result and c(r) so i ar the solar profile extracted from 
BiSON data. The black and red profiles correspond to the high-metallicity GS98- 
SFII and low-metallicity AGSS09-SFII SSMs, respectively. 
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Figure 6: (Color online) Left: A schematic illustration of the MSW crossing 
for a normal hierarchy and small 9 V . The dashed lines - the electron-electron 
and muon-muon diagonal elements of the matrix - intersect at the level- 
crossing density p c . The solid lines are the trajectories of the light and heavy local 
mass eigenstates. A v e produced at high density as ~ vh will, under adiabatic 
propagation, remain ~ vh, exiting the Sun as ~ u^ L . Right: the white "MSW 
triangle" is the region where a level crossing occurs and propagation is adiabatic, 
producing strong u e —¥ conversion. Regions of three possible MSW solutions 
frequently discussed in the 1990s - SMA, LMA, and LOW - are depicted by the 
blue ellipses. 
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Figure 7: (Color online) The three-flavor neutrino oscillation contours resulting 
from the SNO combined analysis: the analysis employs only solar neutrino and 
KamLAND data, but the results are in excellent agreement with the conclu- 
sions from global analyses that include recent reactor- and accelerator-neutrino 



constraints on #13. From SNO Collaboration, B. Aharmim et al. (2012). 
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Figure 8: (Color online) Survival probabilities P Ue for pp, pep, 7 Be, and 8 B 
neutrinos deduced from global solar neutrinos analyses, Borexino, and the SNO 
combined analysis, compared to the MSW prediction, taking into account present 
uncertainties on mixing angles. From SNO Collaboration, B. Aharmim et al. 



(2012), with pep result from Borexino Collaboration, G. Bellini et al. (2012a) 



added. 
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Figure 9: (Color online) Left panel: The expected A^yfa f° r Super-Kamiokande 
as a function of Sm^i and sin 2 612 for two (solid lines) and three (dashed lines, 
013 = 9.2°) flavors. The lines are the contours of constant ^4^N- ^he bhie- 
triangle marks the best-fit parameters from current global analyses. Adapted 



from Blennow, Ohlsson & Snellman (2004). Right panel: Ajy^ values consistent 
with current values of 612 and #13 (red band), plotted as a function of Sm^i, 
compared to SKI+II+III results (la statistical and statistical+systematic errors 
indicated). The vertical bands are the la KamLAND and solar values for Jm 2 ^. 
Adapted from A. Renshaw and M. Smy (private communication). 
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Figure 10: (Color online) Accretion on a young Sun with a thin convective en- 
velope as a function of the mass and metallicity of the accreted material. The 
upper panels show the deduced He and metal content of the pre-solar gas, Yi n j and 
Zj n i, and present-day Zg, determined from the luminosity, radius, and AGSS09 
photospheric Zs/Xs=0.0178 constraints. Truncated trajectories indicate the ab- 
sence of a physical solution. The lower panels compare the resulting helioseismic 
properties of the models to observation (horizontal bands). 
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Figure 11: (Color online) Global x 2 solar neutrino analysis for the late accretion 
NSSM scenario discussed in the text. Contours are shown for 1%, 10% and 68.4% 
confidence, with the best-fit model indicated by the red diamond. The fixed-x 2 
contours corresponding to the AGSS09-SFII and GS98-SFII SSM fits are overlaid, 
showing the very limited parameter space where NSSMs with accretion do better 
than SSMs. 



